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TECHNICAL NOTE 2582 

GENERAL CONSIDERATION OF PROBLEMS IN COMPRESSIBIE 
FLOW USING THE HODOGRAPH METHOD 


By Chieh-Chien Chang 
SUMMARY 


The purpose of the present report is to investigate the hodograph 
method as it is applied in general to the problem of compressible flow. 
First, the hodograph equations are given in various canonical forms 
which are convenient for obtaining solutions in the different flow 
regimes. 

Since the coefficients of the canonical differential equations are 
implicit functions, exact solutions are difficult to find. Consequently 
different approximations are chosen so that some simpler differential 
equations capable of solution can be obtained. For most of the cases, 
fundamental or singular solutions are given or indicated. 

The detailed development is concentrated on Chaplygin 1 s second 
equation. The first-order approximation is well-known as the Tricomi 
equation. The second- and third-order approximations have a rather new 
approach. Both approximations follow the exact gas law closely in the 
neighborhood of the sonic velocity. The solutions are found to be 
Whittaker functions and the associated confluent hypergeometrical 
functions. Both approximations can be applied to the incompressible 
flow so that Chaplygin 1 s procedure of borrowing the boundary conditions 
can be used if necessary. For the third-order approximation, the corre- 
sponding hypothetical gas law is derived and is found to differ very 
little from the exact gas law. The transformation relation between the 
hodograph plane and the physical plane is also given for the various 
solut ions cons idered. 

To make a comparison of the present approximate solution with the 
exact Chaplygin solutions, the flow through an aperture, studied by 
Chaplygin and Lighthill, is reexamined. There is some difference in 
the problem itself, as well as in the method of Chaplygin and Lighthill. 
First, the vessel with straight walls inclined at an arbitrary angle is 
considered rather than that with the wall at right angles. Second, no 
association of the boundary conditions with those for incompressible 
flow is made. The problem is treated directly as a boundary- value 
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problem. The result calculated with .the Whittaker function checks 
well with that obtained by Chaplygin and Lighthill. 


1 - INTRODUCTION 


For the isentropic steady two-dimensional flow of nonviscous com- 
pressible fluids, the pair of differential equations in terms of the 
stream function and the potential 0 are nonlinear, and cannot 
be solved analytically except for a few special cases. Chaplygin 
(reference l) 1 introduced the hodograph method in order to transform 
these nonlinear equations to linear ones, so that the available 
classical mathematical analysis and the principle of superposition can 
be applied. While there is a gain in the linearity of the equations, 
some new difficulties arise in the hodograph method. First of all, 
it is difficult, in general, to transfer the physical boundary condi- 
tions to the hodograph plane. Second, the flow in the hodograph plane 
usually consists of multiple-sheeted Riemann surfaces with a number 
of singularities where the analytic continuation of the series solu- 
tions becomes very complicated. Third, only one kind of a particular 
solution of the stream function - hypergeometrical functions and 
trigonometrical functions - has been obtained so far for the equations 
in the hodograph plane. By superposition, a series solution is 
achieved, but is difficult to apply even if the boundary conditions 
are known. Fourth, the transformation between the physical and hodo- 
graph planes becomes singular when the Jacobian determinate becomes 
zero or infinite. This is likely to happen when the supersonic flow 
is imbedded in a subsonic flow region. The existence of the so-called 
limiting line (reference 7) causes a breakdown of the entire flow. 

In order to obtain some solutions in the case of subsonic flow, 
Chaplygin discovered an ingenious method of obtaining usable solutions 
for compressible flows by comparison with the corresponding series 
solutions of the incompressible-flow patterns as the limiting case in 
the hodograph plane. In the case of no circulation around the closed 
body, the derived series solutions involve a free constant to be chosen 
at pleasure. But fortunately each solution chosen corresponds to a 
reasonable flow pattern in the physical plane. As the free-stream 
Mach number increases, the body shape deviates from the image body in 
the incompressible flow. Thus, in general, the body shape cannot be 
preassigned and the compressible flow about it determined at a certain 
Mach number. This method has been carried further by Tsien and Kuo 
(references 2 and 3 ) > Lighthill (reference L) , and Cherry (reference 5) 
on transonic flow up to the occurrence of the limiting line. Some of 
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their solutions apply to simple closed bodies with circulation. There 
are difficulties in the analytic continuation for an arbitrary flow 
pattern. 

Although Chaplygin and others were pioneers in developing and 
adapting the hodograph method to gas dynamics in the early part of this 
century^ results of the later Russian investigators were not available 
in the English language up to 19^0. Since the second World War, the 
work of Russians has again assumed prominence in the field; namely; 

Frankl (reference 11); Falkovich (reference 34); and Christianovitch 
(reference 42). Along the line of formulation by Tricomi; Frankl’ s 
proof of the existence and the uniqueness of the solution of Chaplygin 1 s 
equation in the transonic regime is interesting. Frankl and Falkovich 
obtained the solution for the channel flow and showed that the stream 
function is triple-valued in the hodograph plane of the axially sym- 
metrical plane at sonic velocity. These results confirm what has been 
achieved by Lighthill (reference 4). Recently; Tomotika and Tarnada 
(reference 33) have formulated some approximate nonlinear hodograph 
equations. A number of interesting particular solutions for the 
channel flow have been obtained. Of course; the solutions are not 
super impos able. Ehlers (reference 10 ) and Carrier; along the line of 
Christianovitch; have obtained the fundamental solutions of the Tricomi 
equation and the corresponding channel flow to the second order of 
approximat ion. 

With the notion of the correspondence of the incompressible flow 
to the compressible; Bergman (reference 6) has developed an integration 
method for calculating the subsonic flow. For the supersonic case an 
extension of the Riemann method is also made by an iteration process 
(reference 4l). Bers and Gelbart (reference 55) have similarly developed 


plex function which satisfies the hodograph equations or the general 
Cauchy- Riemann equations in equivalence. All are important contribu- 
tions to the solution of the differential equation but offer rather 
difficult ways to obtain useful solutions for the flow of compressible 
fluids. 

In the last few years; Guderley (reference 30 ) and Yoshihara 
(reference 36 ) have given a number of papers on transonic flows ; 
particularly with the application of the transonic similarity law which 
was developed independently by Von Karman (reference 29) and Guderley 
(reference 30) . They have achieved very important approximate results 
in transonic flow. 

In the subsonic case; if a linear approximation is used to replace 
the isentropic pressure-density relation; Karman (reference 7) and 


a line-integral operator to construct 
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Tsien (reference l4) , following another approach of Chaplygin, have 
obtained many useful results, particularly the well-known K^rman-Tsien 
formula (reference 7) for pressure corrections. Lin (reference 15 ) hns 
discussed the conditions for obtaining closed contours in the physical 
plane, if circulation exists. Clauser (reference 24) recently applied 
this method to find the body shape with preassigned pressure distribu- 
tion and this is definitely very useful in high-subsonic laminar-flow 
problems. Garrick and Kaplan (reference 20) have another interesting 
approach to pressure-correction formulas using Ringleb’s solution 
(reference 8) of a simple source and vortex. 

The present investigation, as the initial step of the research 
program of this challenging problem, is mainly interested in the 
following three aspects of the problem: 

(a) Besides Chaplygin’s differential equations (reference l), 
can other useful forms of the differential equations be found system- 
atically, particularly the canonical forms in the different flow regimes? 
Of particular interest are the fundamental or singular solutions which 
represent the types of singularities encountered in the hodograph plane. 

Is there any method of constructing such solutions as shown by Picard, 
Hadamard (reference 52), Hilbert, Riemann, and Tricomi (reference 53)? 

(b) If such fundamental solutions are too difficult to construct 
or too complicated to apply to the flow problem, what other reasonable 
approximations, besides the Karman-Tsien (reference 7) approach, can 
be made so that some useful results can be derived in the different 
flow regimes? 

(c) It is well-known that the boundary conditions of the stream 
function + (not the potential I) are well-defined in the hodograph 
for polygonal profiles with or without free streamlines up to transonic 
flow. As a reasonable preliminary approach, can the compressible flow 
around such a given body be found? 

In brief, this paper contains a systematic list of useful forms of 
the differential equations. The canonical forms of the exact differential 
equations are given for subsonic, supersonic, and transonic regimes. 
Unfortunately, one of the coefficients of the equation is an implicit 
function of one independent variable. This makes it impossible to con- 
struct fundamental solutions which would be of value for practical 
application. 

Approximations to the implicit function are chosen in such a way 
that the solutions can be found from classical mathematical analysis. 

In applying the approximations, three objectives are kept in mind: 

First, the differential equation must reduce to the Laplace equation 
as M — > 0 or M > M OT so that Chaplygin' s procedure of utilizing 
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the incompressible flow can be followed. Second, in the transonic 
range, the transonic similarity law of Karman (reference 29) can be 
applied and the simplified boundary conditions can be used. Third, 
the singular solutions and the fundamental solutions in the subsonic 
and transonic regime can be found so that the flow can be determined 
directly if the boundary conditions in the hodograph plane are assigned. 
Or, alternately, the Riemann function in the supersonic regime can be 
found, so that an integral solution can be obtained when the initial 
value or the Cauchy data along a noncharacteristic line are specified. 
Last, the characteristics method can be used. 

In the present paper, a number of approximations for the canonical 
forms in the subsonic flow are given. The zero-order approximation is 
actually the same as Von Karman' s approximation given in reference 7 > 
equation ( 63 ), or reference 13, pages 186 to 188. Both differential 
equations of the first-order approximation can be reduced to the 
Laplace equation in polar coordinates. The singular and fundamental 
solutions have been given. 

The second-order approximations to the differential equations can 
be reduced to Stratton's equation (reference 70) by means of the 
separation of variables. The particular solutions are indicated. 

There are some better approximations which should hold for any 
subsonic Mach number. They are shown in section 3, "Solutions to 
Canonical Forms of Approximate Differential Equations in Subsonic, 
Supersonic, and Transonic Regimes." Similar approximations are obtained 
in the supersonic region. Of course, for each approximation, the 
corresponding pressure-density relation must be determined. 

If the general class of singular solutions in the hodograph plane 
can be obtained, then it should be possible to solve a large number of 
problems by a suitable placement of the singularities in the hodograph 
plane to satisfy the desired boundary condition in much the same way 
that sources, vortices, doublets, and so forth are used in incompres- 
sible flow. 

Some investigations on possible approximations to Chaplygin's 
second equation are made in section 4, "Different Approximations to 
Chaplygin's Differential Equation and Their Solutions." This second 
equation is more convenient to use in flow than the first equation. 

The first-order approximation to this differential equation has been 
shown by Frankl (reference ll) to be Tricomi’s equation. Both the 
second- and third-order approximations are shown to be associated with 
Whittaker' s equation (reference 6 l). The third approximation is suf- 
ficiently good that it should give fairly accurate results in all 
transonic-flow problems. The corresponding hypothetical gas law is 
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shown in section 5» "Hypothetical Gas Law Corresponding to Approximations 
of Chaplygin's Second Equation," and deviates so little from the exact gas 
law in transonic range that the exact gas law may justifiably replace 
it. The transformation relation between the physical and hodograph 
planes is also shown for the third approximation. 

The polygonal body, either closed or open, is known to give simple 
boundary conditions in the hodograph for the stream function but 
not for the potential The particular solutions obtained by the 

principle of separation of variables and the series solutions obtained 
by superposition cannot be applied directly without solving systems of 
an infinite number of simultaneous equations to determine the infinite 
number of coefficients in the series. In this class of bodies with 
one or more convex corners, the solution of the incompressible flow 
cannot be used because the velocity at such a corner is infinite 
according to the theory of incompressible flow. Physically, the flow 
passing such a body is always transonic in character, no matter how 
low the free- stream velocity is. As shown by Guderley (reference 31 ) 
and Busemann (reference 38 ) some shock always occurs at such corners. 

But with the simplified assumption of the transonic similarity law, 
the flow about this group of bodies should be obtainable. In 
section 3 , "Solutions to Canonical Forms of Approximate Differential 
Equations in Subsonic, Supersonic, and Transonic Regimes," a more gen- 
eral simplified equation satisfying the transonic similarity law is 
given and the solutions are shown. 

For an open body built with straight-line elements but with no 
convex corners the problem can be attacked with the Chaplygin technique 
or solved directly as a boundary- value problem. In this case, no free 
constants can be chosen, and the solution is uniquely determined. 

Owing to the complicated nature of the asymptotic behavior of the 
Whittaker function with a very large parameter, and simultaneously with 
a very large value of the independent variable, the future work will 
devote a considerable amount of time to finding the asymptotic solu- 
tions corresponding physically to different flow regimes. This is an 
important step in solving flow problems with the Whittaker function. 

The Whittaker function converges very slowly for a large parameter, and 
is quite similar to the Chaplygin function in this respect. 

In order to demonstrate the method, the flow through an aperture 
is reexamined in order to check with the results of Chaplygin and 
Lighthill. It is found that the third-order approximation agrees very 
well with their work. This is given in the last section of the paper. 

This investigation, conducted at The Johns Hopkins University, was 
sponsored by and conducted with the financial assistance of the National 
Advisory Committee for Aeronautics. 
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2 - CANONICAL AND OTHER FORMS OF DIFFERENTIAL EQUATIONS IN 
SUBSONIC, TRANSONIC, AND SUPERSONIC REGIMES 
2.1 - General Transformation of Differential Equations in 

Hodograph Plane 


The pair of Chaplygin's first-order simultaneous partial differ- 
ential equations for the stream function t and the potential 0 in 
the hodograph plane (equations ( 11 . 10 ) and (ll.ll) in reference 13 ) is 


80 _ Po 1 - M^ 

8q p q 80 


(la) 


80 = Po 8* 

80 p ^ 8q 


(lh) 


where p Q / p and 1 - M 2 are given functions of the independent 

p 

variable q. 

The above system can be transformed to many other forms which are 
more adaptable for analysis if a new independent variable Q is 
introduced to replace q, thus: 


Q = Q(q) 


(2) 


2 

A list of symbols is given in appendix A. 
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In terms of Q, equation (2) yields 


80 dQ = Po 1 - M 2 8f 
8 q dq P q 80 


(3a) 


80 = £0 dQ 8t 

8© p dq 8Q 


(3b) 


Eliminating 0 and then i from the above, a pair of second-order 
partial differential equations is obtained as follows: 


An 2 ’ 2 

8Q q 


- 2 ^ 2,1 


8 2 j + 1 - I^ /dQ\~ 2 8 2 ^ + 8t _d_ 


8© 2 aq 


log c 


/ Po dQ\ 

\ P q dq j 


= 0 


(4a) 


dq 2 q 2 


2 8 2 0 

80 d_ y 

l0 s e 

P 0 1 - M 2 ^ 

QX- 1 

l 

80 2 

8Q dQ 

J* q \d' 

qi 

-- 


■ = 0 (to) 


Of course, if Q is chosen equal to q, the following pair of equations 
is obtained: 


8^ + 1 - 8 2 )|c + 1 + M 2 8i|f _ q 

8q 2 + q 2 80 2 Q ^q 


(4c) 


8 2 0 + 1 - M 2 8 2 0 d /Pq 1 - M 2 \80 _ q 

8q 2 q 2 80 2 d 8 \P q /dq 


(4d) 


which have been used by Ringleb (reference 8). Now make some proper 
choice of the function Q(q), so that either one or both of the above 
equations can be transformed to some simpler form. 
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2.2 - Chaplygin’s Differential Equations - First Form 

In Chaplygin's differential equations of the first form, he 
introduced the algebraic relation between Q and q as 


Q(q) = T 



( 5 ) 


where a 0 * is the stagnation sound velocity. It is obvious that the 

range of r (t has an important physical meaning. It is the ratio 
of the kinetic energy of the gas to the total energy or enthalpy) is 

q 


from 0 to 1. Figure 1 shows the relation of t, 


and Mach num- 


ber M. Equation (4a) transforms to his well-known equation 


8_ 


2t 


(l - t) 


P 8t 


1 - 


di T 


2t( 1 - t) 


P+1 


d 2 j 

8e 2 


= 0 


(6a) 


where 




p. = - + — ' and the relations — = (l - t)^ 3 and 
1 7 - 1 P 0 

are introduced. By means of the principle of the 


separation of variables, the particular solutions of this equation are 
combinations of hypergeometric and trigonometric functions. The dif- 
ferential equation in terms of the potential 0 transforms to a much 
more complicated form 


JL 2t(1 - t) P+1 80 
8t 1 - 3t 

and its solution is difficult to Obtain directly, but it can be obtained 
by integrating the differentials of \|r. Actually, there are four sin- 
gularities of this equation - three regular ones at t = 0, 

t = — = — , and t = 1 and one irregular singularity at r = 00. 

2P + 1 

Therefore its solution is not so simple as the other equation. 


(1 - t )* 3 6 2 0 
2T de 2 


= 0 


(6b) 


10 


NACA TN 2582 


2.3 - Chaplygin's Differential Equations - Second Form 


To eliminate the last term in equation (4a), choose 


Q(q) = cr 


( 7 ) 


where cr is defined by (c CT is a constant) 



(8) 


or 


a 




P dq 



Substituting into equation (4a), the last term drops out and 


(9) 


where 


K d 2 i|r 



(10a) 



In figure 2, K is given as a function of t. Of course, K should 
be expressed in terms of cr. Unfortunately, it is an implicit func- 
tion of a and the differential equation is impossible to solve 
exactly. Figure 3 shows the behavior of K as a function of a for 
the case c^ = -1 and 7 = 1. 4 where the upper limit of the integra- 
tion is q = a*, the sound velocity. It is of interest to observe the 
two asymptotes of the K - a curve: (a) K — 1 as M — 0 or 
cr — ► 00 and (b) K — >• -00 as a — > -0.2513. The function cr can 
be integrated as follows: 
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a 


Os* 


_P_ 


dq 

q 


= -0.2513 + tanh"-*-(l - t) 1 / 2 - (l 



1 - t ( 1 - t ) 2 

+ — 

3 5 


( 12 ) 

where t is given in equation ( 5 ) and y = 1.4 is explicitly intro- 
duced, In figure b, a is given as a function of t. If equation (9) 
is substituted into equation (4b), 



which is rather more complicated. 


(10b) 


2.4 - Another Form of Chaplygin’s Second Differential Equation 
To eliminate the last term in equation (4b), introduce 

Q(q) =V (13) 


where v is defined by (c v is a constant) 

- - 'v 

p 0 1 - M 2 \dq / 


(i4) 


or 


na* 


V = -c. 


p o 1 - M 2 


dq 


( 15 ) 


Substituting in equation (4b), the last term drops out and 
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where K is the same as defined in equation (ll) and c v = -1 is 
chosen. Of course K should he expressed in terms of v, if it were 
not an implicit function of v. No detailed investigation will he 
given in this report. 

Similarly it can he shown that 


It should he noted that equations (10) or (l6) are simpler for only one 
of the pair, at the sacrifice of the other. 


It is well-known that the fundamental solution may he found when 
the differential equation is reduced to the canonical form. Therefore, 
the canonical form is worth while to obtain. If the velocity is sub- 
sonic everywhere in the domain, the canonical forms of the differential 
equations can he found in equation (4) in the subsonic or elliptic 
range. Take 



2.5 - Canonical Forms in Subsonic Flow 


Q(q) = 


(17) 


such that, for M < 1, 



or 



q 


(18a) 


^Actually 
a* can hi 



Here the minus sign is chosen so that 
lit.) After integration, 


can he imposed as an upper limit 


NACA TN 2582 


13 


^a* 


O) = 


(l - M 2 ) 


1/2 


dq 


(18b) 


Uq 


Figure 5 shows m as a function of t, M, and a. In terms of M 
or t , 


u^ S tanh-l(l ^!) l/2 - tanh-Ml - « 2 ) l/2 


0 ) = 


1/2 


= IA-, ' tanh 


Hl T 


^2,(1 - t) 
Equations (4a) and (4b) become 


1/2 


- tanh 


-If 1 - ^1 T 

kl - T , 


.1/2 


(18c) 


d 2 ^ ^ d/. ^lM _ 

5S + S5 + X sK K > - 0 


(19a) 


^ f J-(lcg e K l/2) . 0 

&d 2 ^ 0 2 8m doA ' 


(19b) 


If f = K l/^\|f* a nd 0 = are introduced, the above equations 

can be written as 


d 2 i|c* d 2 ^* 

~ P~ + —k~ + 

dm 8 e 2 


H/rf 1 ir 

(i6\k / ~ 4 k_ 


i|r* = 0 


(19c) 


8 2 0 * a 2 0 * 

dm£ + 8 e 2 


i/K 

’\ 2 1 K" 

_ 1 6\P 

[ / 4 


0 * = 0 


(19d) 


which are the canonical forms of the differential equations in subsonic 
flow and were first pointed out by Bergman (reference 6). 
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2 K 


d(log e K 1 / 2 ) 


-(m - l)mT 2 


(l - m-]_ t ) 3 ^( 1 - t) 1 / 2 


Hi 


( 20 ) 


^ " 1 (l - M 2 ) 372 


*1 " 1 


Figure 6 shows 


d(log e K 1 / 2 ) 


as a function of a). 


Unfortunately the coefficient of the first derivative is an 
implicit function of a) and also is singular at M = 1. This prevents 
the use of Hilbert's and Hadamard's approach (reference 52) to formulate 
the fundamental solution. Bergman (reference 6), Bers (reference 17), 
and Gelbart (reference 27) in this country and Eichler (reference 58 ) 
in Germany have given the integral solution of equation (19) but the 
process is very complicated. The present paper will give some solu- 
tions of the above differential equations ( 19 ) under certain 
approximations . 


2.6 - Canonical Forms in Supersonic Range 


If the velocity q in the flow is entirely supersonic, take 


Q(q) = fi 


( 21 ) 


where Cl is defined hy 



(22a) 


or, if the positive srgn is chosen 



q 


3H 
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After integration, 


a . / aq 


(22b) 


Ja* 


where fi can be expressed in terms of M or t , each in two ways, 


0. = 


Hi 


1/2 


tan -1 ^ 2 ^ — ^ ^ - tan -1 (M 2 - l) 1 / 2 


= Up 


1/2 


tan" 


Hl T 


- 1 


HpU - t) 


1/2 


- tan 


•lPl T 


l\l/2 


1 - T 


(22c) 


or 


= ul/2 


ft = ■- COS' 


Hi 


M + 


K “ X ) 


1/2 


-1 1 

- cos ^ — 
M 


Hi 


1/2 cos-1 


Hid " T )~ 


^1 


1/2 


COS 


-1 


fctVl 


1/2 


(22d) 


Figure 5 also shows fl as a function of t # Actually ft is a dis- 
torted velocity magnitude such that ft = ±0 become two families of 
simple straight characteristics inclined at ±^5°- (Refer to p. 215 of 
reference 13. ) Equation (4) becomes 


^ H -t[[og e (-K) 1 / 2 ] = 0 

8n 2 8e 2 dn l e -J 


(23a) 


8 2 / 820 80 d r n 

— ( - K) j =0 


8ft 8e" 


(23b) 
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In the hodograph plane, the characteristics are the same as those of 
the simple wave equation and are invariant to the behavior of K(ft). 

If t = (-K) and 0 = are introduced, the following 

equations are obtained: 



(23c) 


a> 2 0* a 2 0* 

dfi 2 80 2 



(23d) 


d|log e (-K) 1 ^ 
dft 


K - l)uiT 2 
(n lT - l) 3//2 (l - T) 1 / 2 


^1 M 4 

' 1 ( m 2 - l ) 3 / 2 


(24) 


Figure 7 gives 


a[iog e (-K) 1 / 2 ] 
dft 


as a funccion of 2. The Riemann 


function (reference 52) always exists uniquely and the solution of both 
equations can be constructed theoretically for a region bounded by 
characteristics if the Cauchy data are given along a line that is not 
a characteristic. Owing to the implicit nature of the coefficient 



as a function of ft, the Riemann function and the 


solution will be too complicated to construct. Besides for the super- 
sonic flow the Cauchy data in the transformed hodograph plane are not 
completely known in general. Therefore the Riemann function is not 
of much use in the hodograph method. 


The present paper will also show some solutions to the differential 
equations under certain approximations. It should be noted that ft 
and oo are related by ft = ioo if it is desired to extend the defini- 
tion of oo to the supersonic side. Such an extension seems very 
obscure as yet in its meaning. Fortunately (dft) 2 always occurs in 
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differentiation, and the derivative of the logarithmic term does not 
depend on such an extension. Therefore equations (23a) and ( 23 b) are 
single- valued. 

There is another interesting feature of the canonical form. In 
the hyperbolic region, the characteristics are the same as those of the 
simple wave equation while in the elliptic region the imaginary char- 
acteristics are the same as those of the Laplace equation. They are 
invariant to physical conditions. The physical law influences only 
the first-derivative term. 


and where f(e) is an arbitrary function. (The upper sign is to be 
used for subsonic flow and the lower sign for supersonic flow. ) Then, 


2.7 - Canonical Forms in Transonic Domain 


Now choose Q(q) = e where e is defined by 



(25) 


U 


O 



For M < 1 choose 




Equations (4a) and (4b) can be written as 



1/2 


= 0 


de i2 be 2 de' de 


8 2 0 + f ( 6 , ) _ 80 d 


flog e [Kf(e')] l/2 j 


= 0 


( 27 b) 
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where e' is used for e in the subsonic region. For M > 1 choose 


"\2 M 2 




M - 1 


(25b) 


d 2 * w _«x c> 2 i x ^ d 1 

P p Ac" He" < 

(8e") 2 8e 2 ° e d€ 


log € 


n-Kp 1 / 2 


f<«") 


► = 0 ( 28 a) 


■0? - f(e "> 3 - ^ sp {*• i-^y /2 ) - 0 (*b) 


where e" is used for e in the supersonic region. 

It is interesting to note that the above two sets of equations can 
be combined into a single set if f(e) is an odd function, that is, 



(25c) 


8e 2 


f ( € ) S + 

8e 2 d€ 


-i°g e 


'(-k) 

f(«) 



0 


(29a) 





( 29 b) 


which are valid for both M < 1 and M > 1. It should be noted that 
6 > 0 for M > 1 and e < 0 for M < 1. 


The canonical form of the mixed differential equation is explained 
as follows. If, following Tricomi (reference 52) f(e) is set equal 

to € and -K = (m 2 - l) is introduced, equation (29) reduces to 
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a ft 

de 2 


a e 2 


a 2 t 

e + 

de 2 

a 2 ^ 

e — £- - 

ae 2 


1 a^ 

2 ae 



1 

2 




19 

(30a) 


(30b) 


which are valid for the mixed region. They can be obtained directly 
by setting in equation (23) if M > 1 and -€^ = ^ to 


2 


in equation (19) if M < 1. These may be called the canonical forms 
of the differential equations in the transonic region. The solution of 
the exact differential equation is obviously difficult. Some singular 
solutions of the above equation under approximations can be obtained 


as will be shown later. 


It is well-known that the characteristics are fixed for all linear 
hyperbolic differential equations. But there is an interesting feature 
in the canonical form such as equations (30a) and (30b), because the 
characteristics are fixed and invariant to the function K(e ) or to 
the physical problem from which the differential equation is derived. 
The characteristic equations are 


( e „ 0O ) 2 - | * 3 = ° 

which represent two families of characteristics of cubic parabolas with 
cusps at their points of intersection with the e = 0 axis. 

Figure 8 shows e as a function of o> and fi. Figure 9 shows 

e jlog e as 8 function of e. 

3 - SOLUTIONS TO CANONICAL FORMS OF APPROXIMATE DIFFERENTIAL 
EQUATIONS IN SUBSONIC, SUPERSONIC, AND TRANSONIC REGIMES 


In the last section, the general derivatives were shown of the 
canonical forms of the differential equations in the different regimes. 
At the same time, it is found that the coefficient of the first- 
derivative term in any one of the canonical forms is always an implicit 
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function of one independent variable. The exact solutions are very- 
difficult to obtain. It seems justifiable to seek some differential 
equations with an approximate coefficient for the first-derivative 
term so that solutions are relatively easy to formulate. Of course, 
the validity of the solutions thus obtained must be examined and the 
limits of the range of application must be defined. In this section, 
the solutions of the approximate differential equations will be given 
and their applications will be discussed. 


3.1 - Subsonic Canonical Forms - Equations (19a) 

and (19b) 

In the subsonic case three approximations to the canonical form of 
the differential equation will be given. Figure 10 shows the exact 

curve of K 1 / 2 = ^(l - M 2 ) 1 ^ 2 against o. 

Zero-order approximation of Von Karman .- In reference 7, Von Karman 
suggests taking 


K 1 / 2 = ^(l - M 2 ) l/2 S £2(l - M*, 2 ) 1 / 2 (3D 

K M oo 

where p M and are the density and Mach number at the free stream. 

In figure 10, the horizontal dotted line shows the approximation 

K . It is apparent that the approximation becomes better as 

UJ^UJqo 

— > 0 or M — > 0. With this approximation, the coefficient of the 
first derivative becomes zero in equation (19) and so 


d 2 \|f + 
doo 2 + be 2 


(32a) 


+ = 0 

du£ de 2 


(32b) 


which are Laplace equations. These equations are invariant under 
translation and rotation. Therefore, to can be set equal to 0 at 
sonic velocity a* and so defined as 
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CD 


r (1 . m 2 ) 1 / 2 


Uq 


q f 


dq 1 


(33) 


Of course, Chaplygin 1 s procedure of utilizing the incompressible-flow 
pattern to find the corresponding compressible flow can be followed. 
There is another possible approach. If the types of singularities and 
their locations in the flow of the hodo graph plane are known and if, 
in addition, the boundary conditions of the stream function are known 
under reasonable assumption or under assignment, the solution of the 
flow should be obtained uniquely by superposition of the fundamental 
solutions corresponding to the types and location of the singularities. 
After the solution in the hodograph plane is obtained, the corresponding 
flow in the physical plane should be checked to determine if it is pos- 
sible or not. To serve such a purpose, the fundamental solutions of \|r 
and 0 due to a source at cd q and 0 Q are given as an example. 


</>(&, 6; CD o ,0 o ) = -log e [((I) - 0 * 3 ) 2 + (0 - e o ) 2 ] l/2 ( 3^a) 


CD - 0) 

(%,e 0 )= -tan " 1 — (3te>) 

y - y o 

Of course, with the simple source, there can be built up the potential 
and the stream function of higher-order singularities such as doublets 
and quadrupoles. Actually, it is more convenient to treat the problems 
with functions of complex variables, because \|r and 0 are harmonic 
functions. 

First-order approximation of the canonical form .- It has been known 
for some time that the above approximation is not very good at high- 
subsonic velocity. A higher order of approximation can be made as 
follows. Instead of taking the approximation as given in equation (3l)> 

the approximation can be made accurate to the slope of the curve K^"/^ 
against cd at the free- stream condition, that is, 



a + bcD 


(35a) 
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where 


a = 



aic^N 

CD ] 

^ / (JO^CDqq 


(35b) 


b 


dK^fj 

k d.CD 760=00^ 


(35c) 


if cd = cDqq is chosen as the free-stream condition. This approximation 
is shown in figure 10. Then equation (19) yields: 


b d\|r 

3o> 2 Se 2 a + bco da> 


d 2 0 8f0 + b 80 = 

3co 2 3e 2 a + ^ 

a 

Introducing to = cd + — into the above equations, 

b 

a 2 * a 2 t 1 at 

a«f ae 2 ’ “ 25 ~ 


( 36 a) 


( 36 b) 


(37a) 


8 2 0 8 2 0 1 80 

3o^ 3e 2 ^ ^ 


(37b) 


These are identical to the equations of Stokes and Beltrami for axially 
symmetrical flow. The particular solutions of interest are the well- 
known Bessel and trigonometric functions, and are not given here. 

These differential equations are invariant under translation in 0. 
From Lamb ! s ’’Hydrodynamics , ” the singular solution corresponding to a 

source at cd^ = - ~ or cd q = 0, 6 Q is 


4H 
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0( £D + P e > - P 0 o) = 


H 03 + W + ( e - e °) 2 ] 


1/2 


(38) 


Recently Weinstein (reference 56) has shown that, for the axially 
symmetrical potential, the fundamental solution for a source at 

K + P e o ) is 


0(cd + |,0; cd 0 + J,e Q ) = 


-log e [(cD - %) 2 + (e - e 0 ) 2 J 

* t) 1/2 (» * t ) 1/2 


ni/2 


+ R 


(“ + |,e) 


(39) 


where the function + ^, 0 ^ is regular at the point ^cjd q + 


The corresponding singular solution to the stream function is 


t(o) + £,0; -^,0 O ) = 


0 - 0 , 


[(“ + f f + ( e • 


^o ) 2 ] l/2 


(to) 


The fundamental solution of the stream function is not given here, 
because it is long and involved. The singular solution for more com- 
plicated sources can be built up easily. Therefore, if the boundary 
conditions and the locations and types of sources are given in the 
hodograph plane, the flow in the hodograph plane can be found with the 
classical technique of boundary- value problems. If the boundary con- 
ditions are not known, the series solution of the incompressible flow 
may be considered and Chaplygin 1 s procedure followed. 


Second- order approximations , - Choose the following approximation: 


where 
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Ul/2 ,2 

b = I CD 




dcu dco 2 /co=a> 


_ 1 /d 2 W 2 \ 

2! \ dco 2 / 


This approximation is also shown in figure 10. Then equation (19) 


8 2 t + 

8 2 t 

b + 2 ccd _ 

8cd 2 

8e 2 

a + bo) + ccd^ ^ 


8^ h 


b + 2 ccd 80 _ q 

Sod 2 

3e 2 

a + txD + ca£ ^ 


These equations can be written in another form by introducing 

cd' = — T — (co - — \ (if b 2 > 4ac). 

(b2 - ltac^M 2c y V / 


8 2 ^ + 8^ _ / 1 + 1 \8i = 0 

Sco' 2 8e 2 W - 1 + v^' 


i + ^ + ( 1 + 1 W 

1 2 ^q 2 \a>' - 1 cd' + lydcD 1 


(41c) 

(4ld) 

yields 

(42a) 

(42b) 

(43a) 

(43b) 


The particular solution can be obtained by the principle of separation 
of variables, if it is assumed that 
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,e) = )* 2 (v) (e) = cos ve} V v)(ay ) (44a) 

0(m',e) - 0 1 < v) ( m ’)0 2 (v) (e) - sln v fU (v) («,’) (Mb) 

x cos vej x 

where ^ and satisfy the ordinary equations 


d 2 ^ 

(dof ) 2 


2cd' 


dt x 




(45a) 


- v 2 0, - 0 

(duf ) 2 (of) 2 - 1 1 


(45b) 


which are special cases of Stratton's equation as shown in reference 70. 
The particular series solution can be obtained but it is too long to 
give here. 


Comments on above three approximations. - If the above three approxi- 


mations are plotted in comparison with — llog„ K 

dCD ' e 


: (log e K l/2 ) 


against a) in 


figure 6, they are not satisfactory. The zero-order approximation 


maintains a zero value , although the true value of — (l 

doo \ 


log e K 1 / 2 ) 


at 


o>„ 


may not he zero. The first-order approximation becomes roughly the 


zero- order approximation to -^^log e 


at a) . while the second- 


order approximation becomes roughly the first-order one. It is apparent 
that the above procedure can be repeated to find all the approximations 
to the curve in figure 11, but this kind of approximation is not shown 
here. Some other approximations are given to take care of the asymptotic 


behavior of 


&K* 1/2 ) 


as 0 ) 


0 and 0 ) 


oo. 
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Approximations to singular behavior of -^-(log K 1 ' 2 ).- Figure 11 

dCD '0 * 

shows the curve of — (log„ K"*"' 2 ) against 00 . Some simple calcu- 

dco\ e - ' 

lations will show that lim 00 - 5 ^- (log K^/^J = - i and 

&[**. £2(1 - *)& 


CD - — ► 0 
-> 0 as a) 


00. Therefore, a simple, but 


good, approximation is to assume 


-SLflog K l/2 ) = £L ~ a 

dcoV S e J 2K ~ a) 


(46a) 


in 


The constant a can be chosen in two ways: 

(1) Choose a = — i* if the exact behavior at sonic velocity is 
desired. This is shown in figure 11. 

(2) Choose a such that if the flow velocity i 

the neighborhood of free-stream velocity is desired. If cjd^ is 

infinitely large, it automatically reduces to Karman's approximation. 
This case is also shown in figure 11. 

This approximation is important because the fundamental and sin- 
gular solutions can be obtained similarly to equations ( 38 ) and (39) 

except the one-half power is replaced by the ^ power. 

The next approximation along this line is to assume 


i( l0g e Kl/2 ) 


cd( 1 + bco) 


(46b) 


where a and b are two free constants to be chosen. 
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A few choices can be made: 

(1) Take a = - and b to be such that — = — — — r 

3 2K m .(l + 

the interest is in the sonic and free-stream velocities. 

(2) Take a and b such that —7 - - = ( — ] and 

oo(l + ba3) \2 K>/ cd=cu 

00 

( — ) are satisfied, for example. This case, shown in the 

\2K/ a)= o. 50)^, 

graph (fig. 11) checks well with the exact curve. The differential 
equations are simplified to 


if 


p 

8oo‘ 


2 8© 2 o)(l + bo)) 8 co 


^ - 0 


(47a) 


8 2 0 , 8 2 0 , a 80 _ 0 

So) 2 6e 2 “d + ^ 


(47b) 


Similarly, it can further be assumed that 

/ 

K' ^ a(l + bo) 
2K o)(l + coo) 


( 46c) 


where a, b, and c are free constants. One case is shown in the 

graph (fig. 11 ) with a, b, and c determined from the values of — 

2K 

at 0.5o)oo, O)oo, and 2a ) m . It checks very well with the exact curve. 


3.2 - Supersonic Canonical Forms - Equations (23a) and (23b) 

Zero- order approximation . - Following the approach of Von Karman in 
the subsonic case, there can be chosen in the supersonic case the 
approximation 


(-K ) 1 / 2 = £°(m 2 - l ) 1 / 2 ~ — (Moo 2 - lW 2 


(43a) 
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and 


a = 


la* 


(m 2 - l) l/2 , 

i dq 1 

q 


(48b) 


where q is always supersonic in the domain of the flow pattern. Then, 
the coefficient of the first derivative in equation (23) is zero, and 


ho. 2 he 2 


= 0 


(49a) 


ti. 0 

ho 2 he 2 


(49b) 


which has the solutions 


t(fi,e) = F + (n + e) + G llf (fi - e) 


t v 


. 0(a,e) = F^n + 0 ) + aJn - e) 

which are well-known classical wave functions. The approximation is a 
straight line parallel to the G-axis as shown in figure 12 with ft = ft . 

OO 

First-order approximation . - Along the same approach as the subsonic 
case, it is easy to show 


c^V b' 

8ft 2 he 2 a ' + ^ 


f* = 0 


(50a) 


8 2 0 8 2 0 b' 80 

8ft 2 he 2 a ' + 


= 0 


4- 


(50b) 
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where 


a' = 


(-K) 1 / 2 - a d (- K > l/2 


dft 




b* = 


d(-K) 1 / 2 


dfi 


Jn=n 


Q I 

Introducing Cl = ft + — into the above equations , 

b T 


8 2 t _ 1 8t = Q 

8n 2 8e 2 & 

8 2 0 8 2 0 1 80 

8ft 2 8e 2 ft 8ft 

If the following characteristic coordinates are introduced, 


(51a) 


(511) 


\ — Cl + 0 
jj. = - Cl + 0 

and then, 

a 2 t + 1 /at at\ _ 0 

SX h[± 2(X - |jl ) Sjj. / 

a 2 0 _ 1 /80 _ a0\ = 0 

ax, a^i 2 (x, - n)\a\ a \x) 


(52a) 

( 52 b) 

(53a) 


(53b) 
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which are Euler- Poisson' s equations with p = 0' = - -i- in \|r and i 
in 0. The Riemann function corresponding to \|r can he shown to he 


X 0 ,p 0 ) = 


(Po - ^) 1/2 (P - ^o ) l/2 


(n - X) 


F 

2 1 


( _ ■? ~b 1; a ) (54) 


where X 0 (or ft Q + 0 Q ) 
point of interest, gF^ 


£1 r 

and \i Q correspond to + — and 6 q} the 
is the hypergeometrical function, and 


(x - X 0 )(n - n 0 ) 
(X - - x 0 ) 


It becomes too cumbersome to write the Riemann function in terms 
of the original coordinates ft and 0. This approximation is shown at 
ft = ft m in figure 12. Since the Cauchy data are difficult to assign to 

the potential the corresponding Riemann function becomes useless. There- 
fore, with the above Riemann function $ the stream function at a 
point (ft o ,0) can be found by the Riemann method, if the Cauchy data 

are sufficiently given. In general, even if the Cauchy data are suf- 
ficiently known in the holograph, the Riemann method is rather difficult 
to apply. 

Second- order approx imat ion . - Following the same procedure as before, 
is not difficult to show that 


d 2 j _ 8^{f _ / 1 + 1 \<H 

dfl 2 do 2 \n + 1 ft - 1/82 

^ ^ + 1 _ 1 + 

Zft 2 80 2 Vn + 1 S’ - 


0 


0 


(55a) 


(55b) 
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where 



d 2(-K)l/2 


(57) 


c 




00 


The Riemann function of \|f is under investigation, and is not so 
easy to obtain unless a long series solution is adopted. However, the 
particular solutions are of Stratton's type, but are not given here. 


be made, but the solution to the resulting equation will be more diffi- 
cult to obtain. 

Comments on above approximations .- The above approximations become 

worse as — >• 0 or M ► 1 . Even at very large Mach numbers, they 

are not very optimistic approximations. 


Some better approximation of the coefficient 



can 


3.3 - Approximations to Transonic Canonical Forms - 


Equations (30a) and (30b) 

Take the transonic canonical forms of equations (30a) and (30b) 
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Following Guderley (reference 30) introduce the variable 


= (D 


3 Q 1 
e 3 


£7 

= — (supersonic) (58a) 

fi 2 


■ ( 1 )’ 5 ■ 


or 


( subsonic ) 


(58b) 


It is interesting to note that £ > 0 in the supersonic side and £ < 0 
in the subsonic side of the hodograph plane. The equation £ = 1 corre- 
sponds to the pair of characteristics starting at 0=0 on the sonic 
line. The independent variable 0 can be eliminated in equation (30a) 
and equation (30b). Since the two equations are equivalent in behavior, 
only equation (30a) is treated here. Thus there is obtained 


£(£ - 1 — - - .€$ + — 
*2 3 as de 9 




> + 


d£ 


e 2 d 

n (-K\] I 

18 de 

r ge (— )JJ 


= 0 


Now examine what conditions must be imposed on e — 

de 


(59) 


log. 


(?) 


so that the variables £ and e are separable in the above equation. 
First, assume 


tU,(;) = ti(e)\|r 2 (£) 


( 60 ) 


Equation (59) yields 





9 *1 d6 2 9 *i de 2 de 



(61) 
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/ \ 2 

where K = (-^J (l - M 2 ) is a known function of g . In general, the 

above equation cannot be solved unless the following conditions are 
imposed: 


and 


1 2e_ d ^l 

3 " 3* x de 



a s + Ps 


(62a) 


1 e£ 

9 tl de 2 


1 _€_ ^1 £ _d_ 
9 ^ de 2 de 



a sPs 


(62b) 


where a s and 0 S are individual constants to be chosen later. Equa- 
tion (6l) can then be written as 


d 2 i 2 

5d - t) — + 
dr 


[I - (1 + “s - - »A*2 - 0 (63) 


which is a hypergeometric equation. Its singularities are £ = 0, 1, 

and 00. The general solution about £ = 0 is 


i2(5) = 


? F l( a s>PsJ 2^ 

2 F 1 ^ a s + Ps + 2’ 2 %> 0 £ 


1/2 


► (is | < l) (64a) 


The solution about £ = 1 is 


(1 - 02 ' as ~ Ps 2 f x (I - a s ,| - p B ; J - ag - p B + 1; 1 - 5 ) 


♦2(5) = 


.V 2 , M 2“ as "^s / 

5 (l - 5 ) 2 F l(j L “ a s> 


1 ' 3 S > | ’ a s - 3s + 1. 


* x -0 


(|i-ei<i) (64 d) 
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Following Kuramer (reference 6l), all the 2b solutions about the 
3 singular points £ = 0, 1, and <» can be obtained. The triple- 
valued behavior of for £ = 0 can be shown easily. This 

agrees with the results of Lighthill, Guderley, and Carrier. Under 
the imposed condition, it is not difficult to show that 


tpU) 


Constant 




(65) 


where a s ^ 3 S and c s is a free constant to be chosen later. Also, 
if by definition z s = ^ [Io 8e ft)] , its approximate value zf s can 
be calculated from equations (62a) and (62b): 







c 3(P s -Os) 

c s t 

3(Ps _a s) 

c s € 


(66) 


The exact value of z s can be determined from equation (25a) with 
f(e) = e and K = (^) 2 (l - M 2 ). 

It is necessary now to choose the approximate variable z 0 so as 
to have the same ordinate and slope as z s at e = 0 (or M = l) if 

the flow at sonic speed is of particular interest. It can be shown 
that, at e = 0, 


z 


s 


= 0 



(67a) 


1.05 


(67b) 
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These conditions determine 





^ c S g 

1 - c s e 


(68a) 



(68b) 


where c s Si -0.525. With the approximation to equation (30a), the 
approximate differential equation becomes 


d 2 ^ 2c s dt 

3e 2 de 2 1 ~ c s e ^ 
The singular solutions of this equation are 


(69) 


e 


-3a„ 


(1 


t(e,e) = 



where a s should he chosen to avoid the limiting line in the physical 
flow pattern. It should he noted that equation (69) and its equivalent 
equation of 0 (not shown here) are one order higher in approximation 
than the Tricomi equation. In other words, the Tricomi equation is 
equivalent to taking z Q = 0 for all values of 6. For negative 

values of a s , the singularity of t is at e = — . If a s is posi- 

c s * 

tive, there are two singularities of one at € = 0 and other at 
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4 - DIFFERENT APPROXIMATIONS TO CHAPLYGIN'S SECOND 


DIFFERENTIAL EQUATION AND THEIR SOLUTIONS 


From equations (lOa) and (lOb) 





(71a) 



(71b) 


where K = 



is an implicit function of a. The exact 


solutions are difficult to obtain. Chaplygin (reference l) in his 
researches on subsonic gas jets chose K = 1 so that his equations 
become Laplace equations. As shown in figure 3 this is a reasonable 
approximation if the maximum Mach number in the jet is much less than 
unity, because K = 1 is actually the asymptote for the true K - a 
curve at M = 0. Karman (reference 7) and Tsien (reference 52) extend 
the idea to the linear approximation of the pres sure- volume relation 
corresponding to the free-stream condition, also obtaining a Laplace 
equation and achieving many fruitful results of technical importance. 

In the transonic flow, of course, such approximations cannot be applied. 
Von Karman (reference 29) in 19^7 gave the approximate differential equa- 
tions which are valid when the flow velocity is in the neighborhood of 
sonic velocity and when the body in the flow is thin. He also found the 
transonic similarity law which gives a satisfactory prediction of the 
wave drag for thin bodies. In the hodograph plane, the approximate 
differential equations are of the Tricomi type. The concept of the 
similarity law has been further discussed by Kaplan (reference 39) and 
Guderley (reference 31 ) and extended by Tsien to hypersonic flow. The 
main contribution of the transonic similarity law to the hodograph 
method is: Instead of investigating the flow about a given thin body, 

the flow of a body with the same thickness distribution as the given 
one but of vanishing thickness ratio can be investigated. Because of 
the vanishing thickness ratio, the boundary conditions can be simplified 
in the hodograph plane and are shown to be consistent with the approxi- 
mation applied to the differential equations. Under these conditions, 
the problem becomes a boundary- value problem of the type studied by 


NACA TN 2582 


37 


Tricomi. It is unnecessary to use associated boundary conditions from 
the incompressible flow as done by Chaplygin. After the transonic flow 
of the body with vanishing thickness ratio is obtained the transonic 
similarity law may be applied to find the aerodynamic behavior of the 
given thin body. 

The success of the above approximation encourages the author to 
seek some higher-order approximations. First take a look at the 
K - a curve shown in figure 3. It has two asymptotes: One is K = 1 

and the other is a = -0.2513* It is a monotonic increasing function 
of ct within the range -0.2513 < a < 00 ; at a = 0, K = 0. It is 
understood that 0 > 0 corresponds to M < 1 and 0 < 0 to M > 1. 
The main interest of the investigation lies in transonic flow. There- 
fore, if possible, the approximation should be so chosen as to maintain 
the behavior of the exact K - a curve at and near the sonic velocity 
(a = 0) and at the same time preserve the asymptotic behavior of the 
exact curve as M — *- 0 (a ► ») and M y °° (a — y -0.2513)* 

Just as important, the approximation equation should possess solu- 
tions within the reach of classical mathematical analysis. With these 
few criterions in mind, it can be seen that the Taylor series expansion 
of K about a = 0 is not a favorable choice, although the first 
approximation to be shown below is of this nature. 

I 

^-.1 - First-Order Approximation in Neighborhood of Sonic Speed 

Frankl (reference 11 ) in 19^-5 obtained the Tricomi equation from 
equations (71) by letting K £ = acr and consequently 




( 72 a) 


8 2 0 8 2 0 1 80 

3„ 2 ae 2 0 80 


(72b) 


where a 


'dK\ n ) 0 
— ) £ 9.^2 and 

A a Ja=0 


a 



as given in the 


earlier definition. The comparison of the approximation with the 

exact value of K is shown in figure 3* From the figure the range of 
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validity of this approximation is very narrow, only in the immediate 
neighborhood of the sonic velocity a*. The solutions of these equa- 
tions have been fully investigated by Tricomi in his famous thesis 
(reference 53)- Carrier and Ehlers (reference 57) and also Weinstein 
(reference 40) have further investigated the singular solutions of 
these equations; Tomotika and Tamada (reference 33) also made interesting 
contributions to the singular and fundamental solutions. Recently 
Guderley (reference 31) made an extensive study of the singular solutions 
and showed a very important singular solution corresponding to a family 
of airfoils at sonic speed. Guderley and Yoshihara (reference 36 ) gave 
the flow over a wedge airfoil at Mach number 1. They employ an elegant 
method of attack in solving the problem. 

In the early days, Euler and Darboux obtained the fundamental solu- 
tions of the same equation in the pure elliptic and in the pure hyper- 
bolic domains as pointed out by Tricomi (reference 53)* The importance 
of Tricomi 1 s work is the recognition of the differential equation in the 
mixed domain, that is, partly elliptic and partly hyperbolic, and the 
proof of the existence and uniqueness of the solutions in such a mixed 
domain. Following the footsteps of Tricomi, Frankl (reference 11 ) has 
shown the existence and uniqueness of the solution of Chaplygin 1 s second 
differential equation in the mixed domain, particularly on two problems, 
one being the detached shock wave of a wedge and the other being the 
supersonic jet from an inclined- walled vessel. 


From the comparison of the first-order approximation and the exact 
curve in figure 3 it is immediately apparent that some improved approxi- 
mation should have a wider range of validity than that given by the 
linear approximation. The usual technique of taking higher-order terms 
in the Taylor series would not be particularly helpful because they 
would not improve the asymptotic behavior for large negative values 
of K and for large values of cr. Therefore, choose 


4.2 - Second-Order Approximation 



( 73 ) 


where 



at a = 0 and c can be chosen in 


any one of three ways: 
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( 1 ) 


d4l- 

7 ^ 


■e) 


at a = 0, favorable for transonic range 

00 , favorable 

* -0.2513, 


(2) c = a so that K>> = K — ^ 1 exactly as a 
for transonic and subsonic range 


(3) c = - 1/0. 2513 so that K 2 = K — » -00 exactly as a 
favorable for transonic and supersonic range 


It is interesting to note that, in equation (73), 



finite value, as a — > », and K 2 — » 00 as a — Thus, for any- 
one of the three choices, Itj always has two desirable asymptotes. 


Therefore, the approximation should be fairly good in the subsonic range 
and supersonic range. Of course, if Chaplygin's procedure of using 
boundary conditions similar to those for an incompressible flow is 
followed, the second choice is a favorable one. Cases (l) and (2) are 
shown in the figure. Both seem good in the transonic range. Case (3) 
is not shown. The only known second-order approximation is given by 
Loewner (reference 73). It is also shown in figure 3 for comparison. 


With the above approximation equations 


(71a) and (71b) become 


8^ + aa 8^|f 
^2 + 1 + ca 50 2 


(74a) 


8^0 + aa 8 2 0 + 1 80 

1 + ca ^ 0 2 a( 1 + ca) 8a 


(7^b) 


Equation (74a) can be solved as follows: 

Assume the variables in \|f(a,0) separable and let 

\lr(a,e) = ti^ n \oH 2 ^ n \e) 
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where the positive integer n is introduced to characterize the nature 
of the solution 




cos 

sin 


+ c 4 

(n = 0) 

(nA,0) 

(n ji 0) 


where the positive constant X is so chosen that S ^ n (nX0) is 

cos ' 

periodic for any fixed interval of 0. For n 0, satisfies 

the equation 


d 2 ^ 11 ) 

do 2 


n X aa 
1 + ca 



= 0 


(75) 


Introduce a new independent variable z = 2nX^j ^ (a + 

in the range -i ^ a ^ 00 , the corresponding range for z is 
Equation (75) yields 


Obviously, 
0 ^ z ^ 00 


d2-^(n) ( z ) 

do 2 



-17 


(n) 


(z) 


= 0 


(76) 


which is a particular case of the Whittaker equation with m 2 

. ni /a \^-/ 2 

a reason to be shown later). Here k = — (— ) >0. 

2c \c/ 


i (for 


Before the choice of the solution is made from the known results 
of classical mathematics, a clear understanding of the nature of the 
desired solution when cr — > 00 or when the flow becomes incompressible 


is necessary. As 


and 


n (n) <z) 


must behave 
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"2 


like e from the limiting case of equation ( 76 ). Furthermore 
sin nA,0 

will he just a single term of the convergent infinite 




cos n \6 


(n) , 


series which represents i|r(a,6), to be shown later. The term ( z ) 

should remain finite when a — > 00 . Therefore, \|r-|^ n )(z) ~ e 2 as 
a — > 00 is necessary. 

( n) 

With the desirable asymptotic behavior of \|r^ '(z) at large 

values of a in mind, the only choice is 


i! (n) (z) = wf n) ! (z) 


k '-2 


whose integral representation (see reference 6l) is 



f(k + l)e 2 z 1 
2ni 


'(Of) 



( 77 ) 


where the path of integration is a contour in the complex t-plane 
starting from » just above the real axis, encircling the origin in 
the positive direction (counterclockwise), and returning to the starting 
point just below the real axis. In general, for k > 0 there are two 
branch points, one at t = 0 and the other at t = -z, a point on the 
negative real axis (since z > 0). Thus the path of integration must 
be chosen so that it will not encircle the branch point t = -z. 

The asymptotic expansion for large values of a (not simultaneously 
of large values of n) is 


W n ! (z) 
k,-- 
' 2 


2 z k 


1 + 


5 - ( k - l) 2 ]li - ( k - 1) 2 ] • • • [t - ( k - 5 * l) 2 _ 


s: 7 . 


(78) 
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For the case of large values of n, the asymptotic expansion of 

(z) is much more complicated owing to the fact that both k 



and z, being proportional to a , become large simultaneously. Very 
little information on such expansions is available. 

For numerical calculation, the corresponding series solution is 
desirable. Jeffreys (reference 7l) condenses early developments of 
Whittaker, Goldstein, and Stoneley and puts into compact form the 
relation between the Whittaker function and the confluent hypergeometric 
function. The following discussion closely parallels his work except 
for some changes in notation. For 2m = ±p (\i is a positive integer), 
the Whittaker function can be represented by a combination of Rummer* s 
series when the limiting value is taken. Thus, 



where 



# 



^Fj(a, 7 ; z) being the confluent hypergeometric function. 

( a = m - k + i, 7 = -2m + 1 according to Jeffreys 1 notation. 


) 
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where 


U 


(-m - k + i, -2m + lj zj = -- — ... — ^F^fm - k + p2m + 1; z) h 

V r (-m - k + i) 

r( 2 m)z 2m _ / . 1 _ . \ 

— i TV l F l(- m - k + p- 2 ® + U z ) 

IMm - k + ' 2 


He gives the solution of the limiting case 2m > ±M- • For the present 

case if m is chosen equal to -i (7 = 2 in his notation), there can 

he written 



(79a) 


where 


U,(z) = - 

r(i - k) 


(79b) 


U P (z) = — — ^,(1 - k,2,; z)|log z - F(l) - F(0) + F(-k)l + Uo (79c) 

* r(-k) L -JO 

The function F(£ ) is called the digamma function and is generally 
represented by \|f(£ + l) (reference 7k). 


F(0 = ~|l 0 Se r(c + lj] 


= lim 
n — > 00 



n - 


1 

i + 1 


_1 

i + 2 


1 

i + 3 


i + n 


(30) 


(F(0) = Euler constant (- 7 ) = -0.57722.) 
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For the present case, 


u 3 u>.;r 


2 s (1 - k) (2 - k) . . . (s - k) 


rx r(-k) 


si (s + 1) : 


X 


r - k 


r 


r=l 


+ 1 — 

r s + 1 


(-k)r(l - k) r - k Z_r s + 1 


= y- r(i + s - k) z g 1 s~ 
sTl r( - k)r(l " k ) \^i 


(k I integer) (8l) 

Then, the series solution of equation (76) is 


00 / \ 

t(a,e) = (cjCT + C 2 )(c 30 + Ck) + \ A n W 2 . ( z ) sin (nXG> + On) (82) 

Z . — &?- 75 


where c^, Cg, Cy c^, A^, and a n are constants to be determined 

from the given boundary conditions or from the boundary conditions 
associated with an incompressible flow as in Chaplygin’s method, if 
necessary. Of course, the convergence of the series must be established 
for the particular problem in order to be sure the above representation 
is correct. 


4.3 - Third-Order Approximation 

The K - a relation can be approximated more closely by assuming 
an analytic function of a in the range - ^ < a < » of the form 


ag(l + bq) 
(l + ca) 2 


( 83 ) 
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such that 


(1) 

k 3 = 

o 

II 

W 

at a 

= 0 

(2) 

d£3 _ 

dK = 

a at 

a = 0 

dg 

dg 



(3) 

d% 

. d 2 K 

at a 

= 0 


dg 2 

dg 2 



(b) 

Ko(g) 

also 

possesses two asymptotes , one of which can be 


made to coincide with an exact asymptote of K(a). 


The conditions (3) and (4) determine the constants b and c. The 


ab 


above approximation with Kg — * — = 1 


as 


00 is shown in fig- 


ure 13 • It checks very well with the exact value of K for the range 
from subsonic through transonic up to supersonic regimes. The other 
asymptote, a = -O.2583 as Kg = -00, does not differ greatly from 


a = -0.2513 for the exact value of K. There is another important 
advantage of this choice, because the boundary conditions of incompres- 
sible flows can be borrowed as in the Chaplygin procedure. 


- Introducing equation (83) into equations (71a) and (71b), 


8 2 t ag(l + bg) 8 2 t 
8g 2 (l + eg) 2 80 2 

and 


(84a) 



a.g(l + tig) 8 2 0 
(1 + eg) 2 dd 2 


(2b - c)g + 1 80 

g(l + bg)(l + eg) 8e 


0 


(84b) 


To solve equation (84a ) , assume that the variables are separable and 


/ n \a,e) = ♦ 1 (n) (g)t 2 (n) (e) 
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where the superscript (n) is used to show each solution is related 
to n, a positive integer. 


4r 2 (n) (0) 


sm 

cos 


(nxe) 


C 3 0 + c 4 


(n f 0) 
(n = 0) 


(85) 


and t-^ n )(a) satisfies 


a2 h (n) <■ b») t (n) , 0 


dcp 


(1 + ca) 2 1 


(86) 


With the introduction of a new independent variable z = 2nX -)fffi(i + c <r), 

c^ 

equation (86) can be transformed to the well-known Whittaker equation 






dcr^ 


(87) 


where 


i c - 2b . fa" ^ 

k = -n A - < 0 

2c 2 


(c > 2b) 


(88a) 


m = 


1 2.2 / , s 

1 _ n X a( c - b ) 


^imaginary if n - ~ ^ 


/real If ~ b > < 

\ k l 


(88b) 
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The fundamental system of solutions of equation (86) is 


t! (n) (z) - 


c n W k”m < 2 > * K„«-k!m ( ‘ z) < n + °> 


(89) 


Cj a + C2 


(n = 0) 


where C n and D n are constants to he determined. 

The Whittaker functions can he expressed in terms of the confluent 
hypergeometric functions and Kummer 1 s functions if m / as follows: 


(z) 

k,m 


and 


r(-2m) ..(n) 


r( 2m) (n) 


1 — r\ V- (z) * / i) V-. 

(-m - k + -) r(m + k + -) 


(z) (90a) 


<-> 


r( - 2m) , M<?> (-Z) ♦ , r(£m) , M<?> (-z) 


r(-m + k * |) 


± \ -k,m 


(m - k * |) 


x \ -k,-m 


(|arg(-z)| < | nj 


(90b) 


where the Kummer' s series are 


..(n) / \ 2 

K,±m (z) ' e 2 


-■I ±m+ i r(± 2 m + 1 ) 


00 P ( + 


nZ 


r tm - k + 7? + slz 


)■ 


±m+| 


= e ^ z 


r^±m - k + r(+2m + 1 + s)sl 

jF^+m - k + |,±2m + 1; z^ 


M 


(n) 

-k,±m 


(-z) = 


e^ z 


| ±m+I 

= e z 



) ( - 


z) £ 


Besides, the integral solutions are much more general, and the Whittaker functions can he defined <£ 
uniquely whether or not m = ±i. 


w[ n) (z) = 
k,m 


- rf-m + k + -)t 

2iti \ 2/ 


z_ 

2 


P( 0+) 


(-t) m - k ^(l.|) m ‘ k -ie-‘ d t (91a) 


<-) - 


2nT r (' m ' k + ?) e£ ( ' ZY 


\o+) 


(-t) 


m+k- 


M) 


m - k - 5 e -t 


dt 


(91b) 


In the present case z is real and positive and the contours are so chosen that the second 
branch point t = -z is excluded. For more details consult reference 6l. 


Now the question arises whether both or only one of the two Whittaker functions exists in 
the present solution. This can be determined from considering the incompressible flow as the 

limiting case as a — > » (m — *0). Furthermore, ( z) ^os nk6 iS just the nth tenn ° f 

the convergent series solution. In other words if^z) must be finite as a > Now, 

from equations (91a) and (91b), the asymptotic expansions for large values of a are 


W< n) (z) 
k,m 


z 

‘2 „k 


and 


T (n) 


w -kjm (* z ) ^ (- z )' k J 


[ 'f - ( k - 1) 2 ] t 2 - ( k - 1) 2 ] • • • \f - ( k - s + 1) 2 ] 

^ 

\ mZ - ( k * l)!][ m2 - ( k * I) 2 ] • • • [ m2 - ( k * 5 - g) 2 ] 

s=l 


( 92 a) 


si (-z) 1 


( 92 b) 
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Equation (92a) gives an exponentially damping function of z and equa- 
tion (92b) is an exponentially increasing function of z which is 

divergent as a >00. Therefore, the only choice is D n = 0. Thus, 

there can be written 


iKcr,e) = (c x a + c 2 )(c 30 + cij 


B W.^ (z) sin nX6 
n k,m ' J 

where c^ c g , c 3 , c^, A n , B n , 

given boundary conditions or from the boundary £Onditions of incompres- 
sible flow, if necessary. If A n cos and B n sin are introduced,, 

there can be written 


(v4?m ( z ) cos nW 


n=l 


( 93 ) 


ind X are determined from the 


\|r(a,0) = (c-^a + c 2 ^c 3 0 + c^j + ^ A^^ ( z ) sin ( nX6 + ^ (9 1 *) 


n=l 


It is interesting to note that, when m = +i, equation ( 87 ) reduces to 

equation (76). Thus the second-order approximation is just one par- 
ticular case of the third-order approximation with b = c. Of course, 
when the choice is made that b = c = 0, both cases reduce to the first- 
order approximation, equation (72a). 


Along this line of thought, if there is imposed k = 0 or 2b = c, 
equation (82) becomes 


acr 


K = 




(95) 


(l + ccr)‘ 


The differential equation becomes the normal form of the Bessel equation 


d2 t (n) J 1 , tif 

dz2 V *> z 2 


= 0 


(96) 
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where 


= pl/2 


z = 2 ' nX 


a ) 


(97a) 


m = 


1 n^A-^a 


2c 3 


Its solution can be written as 


(97b) 


*l (n> (2) - M - «c “«Km (n) (|) 

where is the modified Bessel function of the second kind with 

complex order m. The function K m ( n ) ^ j and no others satisfy the 

requirement of boundedness in value as 0 — y <». This can be shown 
by the asymptotic expansion 


^ (n) (f) 


.(nightie' 5 J 

ltfZ 


x , f- [ m2 - (?) 2 ] [ m2 - (l) 2 ] • • • [ m2 - ( s - ?)!] 


r 

s=l 


s: z 


(98) 


There is another interesting feature of this approximation if the 

following conditions are chosen: a = — at 0 — y 0 and c such that 

do 

c = — or K — > 1 as 0 — y Then, 

2 


K = 


aJ ( 1 + | CT ) 

(l * | of 


(99) 
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and 


z = nX 



( 100a ) 


1 kn 2 \ 2 


(100b) 


m 


The advantage of this is that the boundary conditions can still be 
brought from the incompressible flow. 

+ 

The accuracy and usefulness of this third-order approximation is 
discussed in the last section of the paper. 


5 - HYPOTHETICAL GAS LAW CORRESPONDING TO APPROXIMATIONS 
OF CHAPLYGIN'S SECOND EQUATION 

In equation (Tla) , if ( cr) ^ ) is introduced, the principle 

of separation of variables yields: 



( 101 ) 


where 



(102) 



1 - T 

3 



(103) 


da 1 

dq P 0 q 


( 104 ) 
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Since K(o') is an implicit function of 0 , some approximations of 
d-iff^rent orders have been made in the preceding section to obtain 
useful solutions of It is natural to ask what kind of hypo- 

thetical gas laws will correspond to the different approximations. 

It is understood that 0 is just another way of expressing the 
velocity magnitude; consequently it is defined once and for all time 
by the exact gas behavior as shown in equation (103). The variable 0 
is independent of the hypothetical gas law for each approximation 
shown in the preceding section. The problem now is to find the func- 
tional relations of p and p/p D in terms of 0 as the independent 
variable. First of all, by definition, 



(105) 


and the differential form of Bernoulli 1 s equation gives 



If equation (104) is rewritten 


_P_ 

p o 




(106) 


it can easily be shown that 



(107) 


Consequently, if is expressed, in terms of q"^- and. its derivatives 

with respect to cr, equations ( 106 ) and (107) yield 


M 


2 



afq- 1 )! 

2 

_ da 



1 


(108a) 
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On the other hand, equations (102) and (106) give an alternate expres- 


which is the same equation as equation (101), if q -1 is placed for f 
and nX is set equal to 1. If the solutions of equation (101) are 
known, the solutions of q can similarly be obtained explicitly in 
terms of cr. 

If now the solution of equation (109) is substituted into equa- 
tion (106), the density ratio is obtained as a function of a. Dif- 
ferentiating this expression with respect to a as given in 
equation (107), the Mach number is obtained in terms of a. The 
derivatives above involve only differentiation but p(a) has to be 
obtained by the integration 


Of course, if equation (109) can be solved exactly, the relations 
of p/p Q and p with respect to cr will coincide with those already 
obtained from the exact gas law. 

Now, introduce different orders of approximation to K(ct) in order 
to solve equation (109). Then, obtain the corresponding approximate 

solution of q“l in terms of a called q _1 (a). With the term q -1 (a), 


tion (110). Then, the approximations can be compared with the exact 
values, both in terms of cr. 


sion for M 2 in terms of q _1 and its derivatives 


M 2 = 1 



K(a) = 1 - 


TZ (l°ge S" 1 ) K ( ff ) 


(108b) 


Equating the above two equations, a differential equation is 
obtained of q _1 with respect to a, 



(109) 



(no) 


p(a)/p 0 can be obtained from equation (104) and p(cr) from equa- 
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In the preceding section, the best approximation to K(cr) is the 
third approximation 





acr(l + bcr) 
( 1 + ccr) 2 


(HI) 


where a, b, and c are constants already chosen according to certain 
considerations. Substituting into equation (109), there results the 

differential equation of q - -*- corresponding to the approximation 
of q" 1 , 


^(q" 1 ) _ ag(l + bn ) (g-i) = 0 
dcr 2 (l + ca) 2 


( 112 ) 


The proper choice of the solution is 


where 


q _ l = EM (z) 

^ k,m 


EM, (z) 


k,-m 


(113) 



+ ccr) > 0 




k = 


c - 2b pr 

2c 2 ^ "b 


< 0 


(114) 




> 0 


The terms m (z) and M^. _ m (z) are Kummer's confluent hypergeometric 

series (reference 6l) as shown in equations (90a) and (90b) and E 
and F are arbitrary constants to be determined from the following 
boundary conditions: 
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(a) At 

a=0, q = q 

= a* 




(b) At 

n ^(q- 1 ) 

a = 0, 

11 

r — 1 
1 

O ' 1 

Tj 

II 

l 

11 

1/7 + 

T 

da 

da 

\ pq /c=0 

a* \ 2 

/ 

In addition, 

q - -*- must go 

to 00 as 

a — y 00 

(M ~ q 

= 0). 


The above are the correct boundary conditions because it is 
intended to choose the hypothetical gas law to agree with the exact 
one at sonic velocity (a = 0) to the highest possible order of approxi- 
mation and to preserve the asymptotic behavior as a — ► 00 . 


2^ab 


If there is introduced z = z Q = 


at a = 0, 


E = 


a * 2 o w ( 2 o) L 


Vt-m( z o) 


ty + l\0 


£ M k,-m( z °)] 


(115) 


F = 


a * z o w ( z o) 


c - Z °K,m( z o) 


(116) 


where w(z Q ) - Mk, m ( 2 oK,_ m ( 2 o) ' M k,- I n( z o) M ic,m( z o) ls the Wronskian - 
From equation (l04) ; the approximate density ratio is 


1 EM k,m( z > + FM k ; -m( z ) 
Po z o c E<^ m (z) + FM^_ m (z) 


(117) 


which is plotted in figure l4. It differs very little from the exact 
curve. From equations (108a) and (108b), the approximate Mach number is 
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which is shown to compare reasonably well with the exact curve in 
figure 15. 

The values of ^ = (a*)^ and p against a can be calculated 
dp 

from 



[“ 4 ,. < 2 > * 


[ E V . <2) * “*,-.<‘>1 {[“ i ,.** 1 * * ™ k , ..<«)]} 


(119) 


and 


P 



da 


Po 


Since the hypothetical gas law differs so little from the exact 
value, it seems justifiable to use the exact gas law to replace the 
hypothetical gas law if necessary, particularly in the neighborhood 
of the sonic velocity. 

To the author's knowledge, the only available high-order approxi- 
mation is Loewners' approximation (reference 73). It is the basis of 
Carrier and Ehlers' investigation on channel flow (reference 32). For 
comparison, his approximation in gas behavior is given in figures lk 
and 15. His approximation is correct at the sonic velocity to the 
second derivative of the K - a curve. 



For the present second approximation b=c or m = - is 

chosen. When 2m is an integer only one of the series solutions is 
valid, namely, M k ^_ m (z). If the solution M k ^_ m (z) is retained then 

the second independent solution can he obtained by using the limiting 
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value of Mj^_ m (z) as m >- i. It is quite similar to the develop- 

ment of the preceding section. No further details are necessary. 

For the case 2b = c, k = 0. Similar solutions can be obtained 
with slightly different boundary conditions. 

6 - TRANSFORMATION BETWEEN H0D0GRAPH PLANE AND PHYSICAL PLANE 


With the introduction of da = the differential equations 

P q 

for \|r and f from equations (10) and (l6) become 


M , K(a) * 

8a 80 


(121a) 


80 _ _8t 

80 8e 


( 121h ) 


where K(a) = j (l - M^). As shown before, the differential equation 


in \| f is 


+ K (a) = 0 

8e 2 


8a 2 


( 122 ) 


with the exact value of K(a) approximated by 


K 3 (a) = 


aa(l + ba) 


(l + ca)^ 
Introducing ty(a,0) = t^( a )^^) ) 


d 2 t 


da' 


2 - - n 2 X 2 K(a)^-L = 0 


(123) 
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where n is an integer and A, is positive and real; A, is introduced 
here in order to make ^(nA,©) = ^(nA,© + 2 k) periodic. For simplicity 

A, = 1 is chosen in the later treatment. With the introduction of 

z = — jab(l 4- ca) (not to "be confused with Z to be introduced 
c 2 

presently), the general solution of equation (122) has been found to 
be 




00 

(a,e) = A 0 (cr + ff Q )(e + e o ) + A n W^ (z) cos n6 + 


n=l 


00 

X 

n=l 


B n W k"m ( z ) sin ne 


(124) 


where A Q , A , B , o , and © Q are to be evaluated from the given 

boundary conditions for a particular problem at hand, or from the 
solution of the corresponding incompressible flow. The superscript (n) 
is applied to the Whittaker’s function in order to show its relation 
to n. Suppose that all these constants are known and the right-hand 
side of equation (124) is assumed convergent and represents ty( <*,©); 
it remains to transform the results obtained in the hodograph plane to 
the physical plane so that the problem is solved in the physical plane. 
The procedure of carrying out such a transformation is given here. 

First introduce q and ©, the inclination of the velocity 
vector, as 


u = q cos 6 
v = q sin 6 


Po 

where u = — £- = is the local velocity component along the X-axis 

8x p 8y 

drf P n 

and v = = - — — ^ oca ^ ve l° c ity component along the Y-axis. 

Then the total differentials of 0, are 
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d0 = M ox + M dY 
8x 8y 

= q(cos 0 dX + sin 9 dY) (125a) 


di = 


St dX ♦ & dY 

8x 8 y 


= ±L q( -sin 0 dX + cos 0 dY) 


( 125 ^) 


Or, in short, 



(125c) 


where Z = X + iY is a complex variable in the physical plane. With 
a, 9 as independent variables in the physical plane, 


d0 = M da + d0 (126a) 

^ 8a 80 

d\|r = da + — d0 (126b) 

8 a 80 


The derivatives of f can be eliminated entirely if the relations 
given in equations (l21a) and (l21b) are introduced in equation (l26a) 


d0 = K(ct) ^ da - d0 
80 8a 


(127) 


6o 


NACA TN 2582 


Substituting equations (126b) and (127) into equation (125c), 
equation (125c) can be rewritten as the exact differential 


J ry 82 J 8 Z 

dZ = — da + — d0 

8a 80 


,i0 


k(o) |i + do + i I*' 

_ 80 p 8aJ ^80 p 80 


d0 


(128a) 


or 


8Z 

8a 


>i0 


K (a) 

_ 80 p 8a_ 


(128b) 


8z = e££ /_ 8i . Po)t 

80 q \ 8a P 80 


( 128c ) 


The above relations are conformal and uniquely exist as long as the 
Jacobian determinate 


8(x,y) 

8(a,0) 


+ 0 


or 00 


(129) 


In the case of ■ ~ ^ =0 at (X,Y) some singularity may occur in 
8(a,0) 

the physical plane. If {^ 2^1 _ some singularity may occur in the 

8(a,0) 

hodograph plane as shown in reference 1+8. 

Before carrying out the calculation, a few symbols are introduced 
to make the presentation a little clearer. Let 

i|r (0) = (a + a o )(0 + 0 O ) 

tc (n) = ( z ) cos ne - w ^ n) cos n0 (130) 

♦ = w[ n) (z) sin n0 = sin n0 

^ k . ) m 
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Also represent Z = corresponding to and Z = Z^ n ^ to 

\|/ n )(z). For future use, rewrite equation (106) as 



(131a) 


and 


1 dp = Po m 2 
p da p 


(131b) 


d_/£o\ = K 
daypqy q 


(131c) 


dZ = ^ da + d© is an exact differential, 

8a 80 


rv^uauiuu y 

(•• i) 


can be 


To carry out the evaluation of Z in terms of a and 0 is 
rather too long for the space available in this report. The basic 
method of calculation is rather simple. Since in equation (l23a) 

8 Z 
86 

integrated to obtain Z except for an unknown function F(a) (or 
F(0)) to be determined. The unknown function F(a) (or F(0)) can 

be evaluated from the remaining relation — (or — ). Of course, in 

8a V 30/ 

the operation, repeated use has been made of the differential equation 

p 

dfc a) p , . 

- - a ) = 0 or Whittaker's equation. For further details 

da 2 

the following section may be consulted. 


n = 0.- In this case, 


z (o) . x (o) + M 




^(a + a 0 ) - lj - i(0 


O 


+ 4 


(o) 


where Z n ^°^ = X n ^°^ + iY is the constant of integration. 


( 132 ) 
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n = 1 . - For n = 1, the function F^ d )(a) is an integral 
ip(l) 


as 


'(cr) = i 


The term Z can he expressed by 


Po cHr 

- £ w<D 

[pq da 

q 


dq 


(133) 


(1) _ Z ( 1 ) = i®! 


and 


z s ^) - z s = 

b b 0 


ie 2ie 

p w<1) . w u) 

- ±] 

IdW^ 1 ) 

. 22 „<D 

4q 

L da P 

2q 

L da 

p 


Po dw(l) (1)1 

— KW V ; da 



2 J « 

p da 




2i6 

e 

dw<1) . p ° M d) 

10 

>> 


4q 

da p 

' 2q 

_ da 

- — 

p 


(I3^a) 


1 

2 


P Q dW^ 1 ) T „ T (!) 


p da 


- KW'' 


da 


(134b) 


where and Z s ^( d ) are the integration constants. 

n / 1 or 0 .- For the cases n / 1 or 0, the function F^ n ^(a) is 
found to be a constant 


z M - Z c <-) . Hi (n ' l)e 
0 2(n + l)q 


dW 


^ w (n) 


da 


le 


-i(n-l)0 


2(n - l) q 


_ 5£o w ( n ) 


da 


(135a) 
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(n) 


e i(n+l)0 

dW^ n ) 

2(n + l)q 

_ 4 or 

e -i(n-l)0 

dW^ n ) 

2(n - l) q 

— da 






(135b) 


From the earlier definition, 

Z = A 0 Z (o) + A 1 Z C (1) + B 1 Z S (1) + V~ [A n Z c (n) + B n Z s (n ^ 


n=2 


(136a) 


and the constant term is 


00 

Zo - A 0 Z 0 (°> + A^d) + BjZ^d) + (A„Z Co (n) + B n z So <n ^] («*) 

n=2 

Substituting equations (132), (134a), (134b), (135a), and 
(135b) into the above equation, 

po/ \ / Cl (^1 + B l) e210 [dW^ 1 ^ 

P ( ° + °o) - 1 - it 8 - So)] ^ A ^ 1^- * 


z - Zo = 


Ao el6 Po 


£0 (1)1 _ ( A 1 + iE l) e ( dW^ 1 ) Po w (l)"j 

P J 2q [_ do p J 


iAl + B 1 | 1 


nwd) 

p da 


da + 


00 

n 

n=2 


~(iA n + B n )ei(n+1)0 

™ {n) + *>o „(n) 

2(n + l)q 

da p 



(-iA n + B n ) e -1 ( n-1 )® 

dW (n) n Po 


2(n - l) q 

da P 


W 


(137a) 
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In more detail, the expression is 


Z - z 0 - ^ 


i0 


+ ®o) - 1 - 1 ( e + ®o) 


(,) 

p k ; m 


(A], + iBi) 0 


2q 


,(D 


da 


(iAi + Bi)e 2ie 

dw kU < z > 

4q 

da 

■ - P ° w' 1 ' (z) 

p k,m 

+ 


n 


iA^ + 


Po dW k^m ( 2 ) 
p da 


K(a)w[^ (z) 


da + 


n=2 


(^n + B n) e 


i(n+l)0 


2(n + l) q 


dwj 1 ^ (z) nD , , 

- - ^ — * — i n l w 

da p k ; m 


(-iA n + 

<*> 

2(n - l) q 

da 


nP ° (z) 

— W k,m ( z .> 


(13Tb) 


2n 


where z = — \Jah( 1 + ca) 


The above result is also true for the cases of the second approxi- 
mation b = c ^m = - or of the approximation b = 2c (k = 0). 

7 - FLOW OF COMPRESSIBLE FLUID THROUGH AH APERTURE OF A TWO- 
DIMENSIONAL INCLINED-WALLED, STRAIGHT-EDGE D NOZZLE 


To apply a critical test to the present investigation it seems 
reasonable to compare the present approximation with a well-known flow 
pattern which has been studied by early explorers using the exact method. 
For instance, Chaplygin gave an application of his investigation to the 
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efflux of a gas from an infinite vessel. Later Lighthill, following his 
approach, repeated the example. Both considered the particular case 
with the walls normal to the jet and both brought the boundary condi- 
tions from the incompressible flow and treated the maximum jet velocity 
up to the sound velocity. 

The present application gives a more general type of such flow - 
flow through an aperture of an inclined- walled, straight- edged nozzle. 
Besides, the problem is treated directly as a boundary- value problem. 


The question of the maximum velocity in the jet depends on the 
value of the ratio of the pressure p surrounding the jet to the stag- 

7 

P fy + l\ 7~1 

nation pressure p Q in the vessel. As long as — ^ ( — - — ) , the 

Po \ ' 

maximum velocity can never exceed the sound velocity. The boundary 
value of the stream function t is clearly defined; therefore it is a 

7 

P h + i\ 7"-*- 

direct boundary- value problem. For the case ^ — J the 

velocity of the jet will be supersonic and this problem is as yet 
unsolved. 

The approximate differential equation used is (given in equa- 
tion (84a)) 


8% aa(l -t- bg) 

(l + c a ) 2 dO 2 

and its general solution has been found to be 


( 138 ) 


\lr(cr,0) = (f^cr + C2)( c 3 e + c^.) + ^ (A n cos nX0 + B n sin nXfi) (z) 


n=l 


= (c^ + c 2 )(c 3 e + c^) + ^ A n W ^ (z) sin (nX6 + Oq) 


n=l 


(139) 
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where c^, c 2 , c^, 

3 rid. Aj^ cos cc^ = 


c A n , B , and X are constants; A n sin = A n 
are obvious; z = ^ab(l + ca) ^ 0, 


k - - DX(c - 2b) J° < 0, and 
2c 2 I* 

tive integer. It represents 
converges. 


in = 


I I (nX,) 2 a(c - b) ^ 

— - — ; and n is a posi- 


the required solution if the series 


The constants c^, c^, 

from the boundary conditions 


c~, Ci . A , a , and A, 
3 4' n' n' 

which are shown as follows. 


are determined 


0 o arbitrary.- The flow in the physical plane is shown in fig- 
ure l6(a.) and in the hodograph planes in figures 16 (b) and 16(c). Since 
it is known that the discharge Q from the aperture is finite, the 
stream function ^ is bounded. Moreover, the flow is symmetrical with 
respect to the center line. If the rate of total discharge is introduced 
as Q, then for inside the vessel there can be written 


i at e = e Q 

t = 0 at 0 = 0 ► 

i|r = § at 6 = -e 0 




Thus it is obvious that is an odd function of 6. Corresponding to 
a constant p^, q and cr on the outer surface of the jet are constant 

from consideration of Bernoulli 1 s equation. Thus \|r = - at 


\ 


0 = o Q = Constant 




0 < 6 <? e Q 


and 


(l4l) 
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First of all, since \|r = 0 at 0=0, 


“n = 0 
C4 = 0 


(lte) 


Thus equation (139) becomes 


i(0,e) = (c-^ + c 2 ) c^e + A n W k^m ( z ) sin nX9 (1^3) 

n=l 

which is an odd function of 0 as required. Next, the stream function 
is defined in the interval -0 Q ^ 0 < 0 Q or. the period 2 0 Q . Thus 

(when n is an integer) sin n\(d + 2 0 O ) = sin nA.0 or 2nX,0 o = 2n«. 

x = i- (ikh) 

9 o 

Substituting equation (l44) into equation (l43), 

00 

f(a,0) = (c 1 a + + S~ A n W^. n ^ (z) sin ^ 0 (l^5) 

fe ' 0 

Now, since is bounded everywhere, particularly when a — * 00 
or M — »• 0, c^ must be equal to 0. Thus equation (l45) yields 

00 

«<T,e) - A o e * A n w k?m < z > sln f 1 6 
n=l 0 


( 146 ) 
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At 9 - 9 q , \|r - - 


Thus 



•^o^o 



(1^7) 


*(,,8) - - ^ A n w[^ (t) sin 0 


( 143 ) 


n=l 


Finally the constants A n can be determined from the boundary condi- 
0 

tion ¥ = - — at a = a n or z = z~ = 


2nn 

V 2 


^ab^l + ca Q j for 0 <9 ^ 9 Q 


* 


' ‘1* 'If * f v4?i (*o) (149) 


° n=l 


or 


Q 0 \ \ . T ,(n) / \ . njtO 

- 51 1 - Tj = 2_ A n«k,m (- 0 ) »ln 

N 0 ' n=l 0 


Now the left-hand side of the equation is a known function 0/6 q and the 
right-hand side is the Fourier series of fifl - ±A i- 

2\ ej 

verges. It is easy to determine the coefficient as 


if the series con- 


A n w k”m W " - 


"II 

-1 


|(l - £) sin (n, l) d(X) 


mt 


( 150 ) 
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or 


A = _ S 1 

n nn t T ( n ) 


W- 


k,m 


( z o) 


Substituting into equation (l48), 


t(a,e) 



w k?il ( z > 

Tf (n) 

nW k,m z 0 


sin 



(151) 


First, it is necessary to show that equation (l5l) is valid for an 
incompressible flow which can be considered as the limiting case of 
a* >oo # From the definition, 



there results, for incompressible flow, 



lim log e 

a* — >00 


a* 

q 


(152) 


and, from the definition of z for very large values of cr, 


^ since ^ = lj 


2nX / , \ 

z = ( 1 + ccx) 

c 


= 2nXff [l + o(ff-lj] 


From equation ( 87 ) the asymptotic expansion of (z) for large 

values of a is 
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*3 

Then the ratio of the two Whittaker functions in equation ( 151 ) at very 
large values of a is 


= e" nXcT (2nX.a) k 1 + 0^ 


4:1 m 


w. 


(n) 

k,m 


( z o) 


= e' nk(tJ_cr o ) ^ k ^ 


0 


$ 


(153) 


Now, for an incompressible flow with equation (152), 



= -log e 


_q_ 

^o 


and 


lim — = 


a*=a r 


a 

C— 


lim 

a*=a 0 *- 


log e a* - log e q 
!°ge a o* - lo Se % 


= 1 
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Introducing the above results into equation (l53)> 



i 


Thus the stream function for the incompressible flow is 



Q JL + Q 

2 0 O rt 



( 15*0 


which reduces to Chaplygin* s case if 6 0 is set equal to tt/2. 
The convergence of the series is easily proved, because 


decreasing function of z in the subsonic range (see appendix B). 


term of a convergent series, and each corresponding term is smaller. 
Therefore, equation (151) is convergent. 

The main interest of this problem has been the minimum width of 
the jet which occurs at X = 00. Therefore the solution must be trans- 
formed back to the physical plane. 

With the differential relation from equations (la) and (lb), the 
following equations can be derived: 



^ 1 for z ^ z Q and W-j^^ (z) > 0 is a monotonically 



M = K * 

dcr 80 


> 


(155) 


dj 

h6 ha 
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and, rewriting equation (128a), 


dZ = dX + i dY 



Differentiating equation (l5l) with respect to a and 9 and 

bering z a = ^ = ^£2 \|ab = ^22 /since = l), 

a da 0 o c * 0 O V c I’ 

St . dz 

8a dz da 



(156) 

remem- 


(157a) 


(157b) 
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Substituting the above equations into equation (156), 


>i0 


-dZ = 




1 + y 

2 Z, .,(n) 


{n l M 


COS 


nit 0 


-1 («o) 


g p 0 g 

P®n 


W V ( Z ) 

v m . n«6 

sm 

®0 


”=1 (M 


dcr + 


>i0 


2Q y“ k,m 

' e o L y(n) 


W, ^ U) 

v m . njt0 

— sin + 


~ r W v 
n=l \ 


% ( Z o) 


0, 


iP 0 Q 


P0 C 


l ^ — w| n) (z) 

i + - ,m 


n=l ( z ) 

k,m \ °/ 


cos 


mt0 


d0 


(158) 


which is the conformal relation of any elementary length in the hodo- 
graph plane to the corresponding element in the physical plane as long 

^ ( X Y ) 

as the Jacobian determinate — - — - — - ^ 0. This condition is automatically 

8(a,0) 

fulfilled in the subsonic region up to sonic velocity as has been shown 
by Tsien (reference 45) and Craggs (reference i*8). Since dZ is an 
exact differential, there can be introduced. 


_bZ 

bo 


e 


ie 


,QK 

b 


w, (n) 

1 , \ W k,m 

2 ^ ,,(n) 


(z) 


n=l W. 


(n) (z ) 

k,m V °/ 


cos 


nit 6 


00 T 


12 PpQ 

P®o b_Y W,( n ^ ( z 0 ) 


ySi 

W^ n - 
n ~l k,m 


(n) (z) 

m ( ^ sinM 

0^ 


(159a) 
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and. 


az _ e ig 
a© q 


V- 2 - it'U) , \ 

2Q \ ^k, m ( z ) _ _ nit© 


L e ° n=I H (■»> 


sm + 

6 n 


iP 0 Q 


P Q c 




»k"m ( z > 


I “- 1 "k?i ( 2 o) 


COS 


njt© 


(159b) 


az 


The value of Z can he determined by integrating — , leaving a func- 

O0 

tion of a to be determined, and, consequently, 


-Z = 


Qe ie Y ***>' 
G.n ^ — r* (n 


( 2 > <*) 


°° q ^ a (, o) 


not nit n . . nit n 

r cos r 8 - 1 sm r e 
u o u o y o 


m) - 1 


QPo e10 

V^ w k?m ( 2 > 

. nit 

nit n 

sm — 6 - cos 

0 O 

£* 0 
e o 

, QPo eie 

0 O P1 

W S ( z o) 


(€ - ■ 


2&bM 


+ F(ff) 

( 160 ) 

where F(cr) is to be determined. This can be done by differentiating 
equation (160) 


az _ Qe 10 Y"~ 
aa © 0 q Z__ 


cos 


mt© 


n=l /nn 



- 1 


w k"i ( 2 o) 




1 lz > - r 'i% < 2 >] 


© - 


sin jlA 0 

"n 


° 9 n ' 1 - 3 m m> (Zo) 

~(2*f - ll— w' (n 

\ 0 o/ J P k ' ra 


/(n) 

‘(rkjm 


-2z rt Wu'“ / (z) 


^ (z) + £21 Kw[ n ^ (z) 

© 0 k,m 


+ K + F’ (a) 

2 0 o 9 


(161) 
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, 2nn 

where z a = - — , 

equations (l59a) 


and equations (l31a) and (l31b) are used, 
and ( l6l ) , 


Equating 



(z) - 


2 2 
nrn 


KW 


(n) 

k,m 


(z) 


W, 


(n) 

k,m 


(* 0 ) 



mr 


do. 



sin 


mr 9 r— 


W- 


(n) 

k,m 


( z o) 


J k,m 


(z) - 


n 2 * 2 (n) 

6cT 


(z) 


= 0 


(162) 


The expressions in the brackets of the right-hand side of equation (162) 
are identically equal to the differential equation, equation (78) if 


2 


V 


2 2 
n it 



and consequently are zero as shown. 


Thus, choose 


F(a) = Constant = Z Q 


(163) 


Now equation (l60) can be rewritten in a more suitable form 


Z = 


-PpQ 

2p S 0 q 


? ie + 


OO 

E 

n=l 


(n) 


(n) 


2 ^ W k,m ( z ) + W k,m ( z ) 




w| n) (z ) 
k,m \ 0 / 


2 p w'^ ( 7 .) (z) -if— -l)' 

2 w k,m U) - w k,m ( z > e \0 O ) 


4*1 ( z o) 


njt _ 1 


njt 


+ 1 


(l64a) 
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Specifically, 


X = 


-PqQ 

2 P0 o q 


<cos 0 + 


00 

L 

n=l 


2 g- w! (n) (z) + w[ n) (z) cos (2i + l)0 
P Q k,m k,m \ e o J 


w[ n) (z ) 
k,m \ 0/ 


2 ife M k ( ,m < z > - W k?i < z > cos fe- l ) 9 


w[ n) (z ) 
k,m \ 0 / 


rut 1 
0o " 


rut 

— +1 

e o 


(161a) 


Y = 


-PpQ 

2 P e o<l 


<sin 0 + 


00 

Z 

n=l 


S ^ U Km < z > + w k?m < z > sl " (^ + l)« 


T U) 


W k"m ( z °) 


21 + 1 

0 O 


2 £ v k,m < z > - 4“i < z > sin (t; - i) 6 


r( n ) 


( z o) 


rut ^ 
0 O ' 


(164c) 


The convergence of the above two series must he established. First, at 
a = a Q or z = z Q , equation (l64b) does not converge for 0=0, because 
the velocity vector of q = q^ at 0 = 0 is located at X = °». 

At very large values of n, 
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( n ) 

With the above notations the asymptotic expansion of m (z) for 
large values of n can be introduced as follows: 


W- 


(n) 

k,m 


(z) = 


■ h '©'(> *§-■) 


,!/6 




ze 


4k- 1 h 


i + + °(si 


(166) 

at sonic velocity (cr = 0). For the subsonic velocity (a > 0), a similar 
expression should be given; the case a = 0 is more critical. Further- 
more from reference 6l, page 352, example 3, 


W,' (n) (z) 


k,m 


v i n i < z > 

k y hl 


k 

z 


1 

2 


E 2 - (* - 


I) 2 ] w L n U <*> 


(«) 


(167) 


Now it is necessary to show that, for large values of n, the value of 

w ’( n ) ( z ) 

k?gI does not depend on n to the first order. Substituting 

W* n) (z) 
k,m 

equation (l66) into equation (167), 


(j] 


1 p - (k - I) 2 ] r(k - i - m) _ s 


= k - - - 

2 


■(k + I - m) 


E + °(=3 


- k - i + (m + k)e' 1|m |r + o(i| 


(168) 


which is a constant independent of l/n to the first order for large 
values of n. Therefore the coefficients involving Whittaker's function 
for large values of n can be written as 


(n) 


2p ^k,m 


( z o) 


P ° W (n) ( z \ 
W k,m V Z ° ) 


| + ( m + k)e~^“J + 1 


+ 1 


(l69a) 
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w' (n) 

2p w k,m 


( z o) 


v l°l ( z °) 


i~ — 


k - | + (m + k)e”^ m 


5 ]- 


(169b) 


It It 

both of which are bounded. For the case Q n = — . — , . 

0 2 3 

an integer) the series in equation (l64c) is convergent, 
shown at the surface of the jet. The equation 


it / . 

. , — (s is 
s 

This can be 


Y = - 


sp 0 Q 




2itpq 

00 00 


sin 6 4- 


2 

n=l 


2 w' / 

2 P^ W k;m ( 

z °) 

1 "I 

sin (sn + 1)0 

J 

N 

O 

) 

sn + 1 


p*, • (n) 

2 P 0 Wk ^ m ( Z °) 

- i 

sin (sn - 1)0 

M (z„) 


sn - 1 

k,m V o ) 


-J - 


(170) 


is dominated by 


S P 0 Q 

2 *PcAo 


1 + E 


^ sin n|0| v sin n|0l\ 

2 + F 2 \> where E 


sin n 1 0 I 


n=l 


and F are positive numbers such that E ^ 
'(n) 


n=l 

; W k ( ,i > ( Z °) 

00 <:i (*») 


+ 1 


and 


F 2- 


'W l \ 

2 ^00 k ,m ' Z o) 


p o w (n) 


- 1 


K) 


for all large values of n. 


series 


Actually fewer terms occur in equation (170) than in the dominate 

r sin n| G I ^ I © I , * 

— = 2 ~ “ — • Therefore equations (170) and conse- 


nt 


quently equation (l64c) are convergent for a = 0. Similarly, the 
convergence for a > 0 can be shown. 
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= 75 - - This is the case which has been treated by Chaplygin and 

Lighthill. The condition 6 = 9 Q occurs at the mouth of the aperture. 
Let h/2 be the half width of the mouth. Then, 


Y = — 


p g Q 


"Pcolco 


n=l 


/ 1 , 1 ) 

1 ( nj 

,PA,m 1 

( z o) 

\2n + 1 2n - ij 

' tr (n) / \ 


J W, 

o k,m 


(-i) n y- (-i) n 

Z 2n + 1 Z 2n - 1 


n=l 


n=l 


p q Q 


^co^oo 


1 + 8 


00 

Z 

n=l 


. 1 ( n ) / \ 
(-l) n n Poo k,m ( z o) 

- 1 (- 0 ) 


(-I)”'- 1 
2n ' - 1 


OO 

r 


(-D 


n-l 


rr 2n - 1 


n=l 


(17D 


If, at X = 00 , the width of the jet is h^, then, by definition, 

P q 

Q = -2^-22 h 


(172) 


Substituting this value into equation (171), 


_h_ _ 2 

h co * 


00 


’(n) 


-l) n n W k,m ( 2 o) f (.i)"- 1 


S3. ^ - 1 P ° ( z o) fo 


f- (-1) 11 - 1 . f- (-I )"' 1 

fes 2n - 1 fcl 2,1 - 1 


1 


( 173 ) 
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Making use of the result of the well-known Leibniz series 
00 



equation (173) can te written 


h_ = 2_ + y~ (-D n n ( z o) + 2/n _ x \ + 2/« \ 

* 11 p o SI 4n 2 - 1 W [ n) (z] 

"k,m \^oj 


1 ♦ ±§s=y 
* p °fci 


( -l) n n W k^ ( 2 o) 

4n2 ■ 1 w ki ( z o) 


(m) 


Just as a check for the case of incompressible flow, the limiting case 
of a* — > co or a — > 00 may be considered. 

First, considering equation (167) as 0 becomes very large and n 
is finite. 


w: ( » } ( z ) 


k,m 


»$, (*) 



nr 


(> - If e-t ^ 



(175) 


The above equation tends to -i when 0 — * 00 , owing to the fact that 

it and m both become zero as 0 — >00 (z — >00). Substituting the 
above equation into equation (17*0> 
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jl . 1 + sf' (-1)°'^ 

h °° 11 n^l ^ n2 “ 1 


= 1 + 


00 

m 


n=l 


(-1) 11 " 1 2 

- — — + _ 
2n - 1 jt 



n'-l 


(-D 

2n' - 1 


+ 1 



(176a) 


Or, inversely. 



— - — = 0.6110 

2 + jt 


(176b) 


which checks with the well-known result of Kirchoff. 

For the case a Q = 0 ( = l) the width ratio h^/h has been 

calculated from equation (17^-). It is found to be 0.7^6 by summing 
terms up to n = 10. It is expected that the error is of the 
order -0.002. Consequently, it checks reasonably well with the Lighthill 
result (h^jh = O.jkkjj which is based on the exact Chaplygin function. 


The Johns Hopkins University 

Baltimore, Md. , April 3, 1951 
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APPENDIX A 
SYMBOLS 


a* 


V 

a. y b y C y d 

A,B,C,D 
pF-^(a,p; z) 


velocity of sound 
stagnation velocity of sound 

constants 

constants 

confluent hypergeometric function 


gF (a,b,c; z) hypergeometric function 


h 

h 


width of jet 

width of jet at X = » 


K ■ ( t ) 2 ( 1 - ” 2 ) 


K-j^ = a a 


Kp = — — 

^ 1 + ccr 

~ _ aa(l + bcr) 

3 = (i / 

if . £ 

z 

M local Mach number 

Moo free- stream Mach number 
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Kunnner function 


m 



z 

n 

P 

Po 

Poo 

Q 

Q 

q 

u 

V 


positive integer 
pressure 

stagnation pressure, pressure in vessel 

free-stream pressure, pressure surrounding jet 

arbitrary function of q (Q(q)) 

rate of total discharge 

local velocity 

X component of velocity 

Y component of velocity 

Whittaker function 



X,Y 


coordinates in physical plane 


Z 


complex variable (X + iY) 




c 


8k 
a, (3 



y ~ ^p/^v 

e = Q(q) 



M. 

7 + 

Vi = 

7 - 

V 


me a tn 2582 

constants 


(defined by * (f^f = - 


2 = q2 

3 = fi2 

2 = e£ 

3 = cn 2 


(supersonic) 


(subsonic) 


inclination of velocity vector 

coordinate of source location, hodograph plane 

slope of nozzle vail 

characteristic coordinate (IT + 0) 


coordinates of source location in hodograph plane 
characteristic coordinate (-JT + 0) 


On* 


Po 1 - M 2 

p q 


dq 


V 


positive constant dependent on boundary conditions ( nX ) 
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P 

Po 


local density- 
stagnation density 
free-stream density 



Po dq 
P q 


7 ~ 




0 

* 

fl = 


potential function 
stream function 


- 1) 1/2 


dq (supersonic) 


— a' 

n = q + — 

b* 


Pa* 


co = 


(1 - m 2 ) 1//2 




CD^ 


dq (subsonic) 


coordinate of source location in hodograph plane 


CD 


cd in free stream 
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APPENDIX B 


PROPERTIES OF WHITTAKER FUNCTIONS 


Here are. presented a few important properties of the Whittaker 
function which are useful in the present investigation. Some of the 
properties are given hy Sharma in reference 59 • He shows that: 

Theorem I: 

The functions W k ^ m (z) and W k ^ m+1 (z) cannot have 
a common zero (root). All roots of W k ^ m (z) are simple; 

W k ^ m (z) and (z) cannot have a common root. 

Between any two consecutive zeros of W k m (z) lies one 
and only one zero of W k _^ m (z). 

The proof is rather simple. See the reference. 

Next, consider the Whittaker equation: 


where n is a positive integer and X is a positive constant. 



(Bl) 


where W = W 



In the present case, 



(B2) 


* = nx(| - i) < 0 


(B3) 
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nr 


nr 


* I • - x ) < ° f ° r n2x - 2 > 

*» A2 < 


4(a - c) 


4(a - c) 


(B4) 


P 1 

In this case mr - — < 0; the powerful Wintner theorem (refer- 
ence 7^-) which requires k < 0 and m^ - i > 0 to make 

4 

/ z \ 

(-l) n ^ 0 (0 ^ z ^ 00 ) cannot be applied. Actually, 


dz 


(z) may oscillate in the range - i < o’ < 0. But the monotonic 
positive nature of (z) for a ^ 0 can be proved as follows: 

Theorem II: 

For nA. > 0 and 0 < a < » (l ^ M > 0), then 
w k)m <*> >°- 
The proof is as follows: 

If equation (Bl) is multiplied by = 4n 2 X 2 ^^, ^^ ere can 

d 2 W n 2 A. 2 aa(l + bn) 


written 


da 


(1 + ca) 2 


W = 0 


or, more properly. 


d 2 W = n 2 A 2 aa(l + ba) 
da 2 (1 + ca) 2 


W 


(B5) 
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By the given conditions, as long as a ^ 0, 


2 

n 


2 

\ aa(l + bcr) 
(l + ccr) 2 


^ 0 


(b6) 


d^W 


(n) 


Therefore — ^ and W have the same sign. But W = (z) is 


da 


O V 

asymptotic to e ^ z = e 
values of a and tends to 0 as a 


./ n ^ + ^nx(l + aj] nX (-t^) 


> 0 for large 


r(n) 


0. Hence W, VXA/ (z) curves 
k,m 


downward as cr 


oo and so must continue to do so while W* 


-1 d 2 W 
da 2 


> 0. 


This proves the theorem. 


There is another interesting feature of the Whittaker equation. 

p 

As long as m^ is real and the function representing the boundary con- 

( n) 

dition is real, W^ m (z) is always real, no matter whether m is 

imaginary or not. For the present investigation, m is Imaginary, hut 

( n) 

^k,m ( z ) * s rea l (0 < z < oo). This question has puzzled the author 

for some time. 
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Figure 1.- Plot of 
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' (l - t ) ; 


and Mach number M as a function of t. y = 1.4. 
^ l/2 

tf = (jttt) * M = 1 > ^7 = °* 9192 at T = !/6; 

M = 1/2, = 1.200 at T = 0.286; M = 1/5 , -1. = 1.583 at t = 0.5; 
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Figure 2.- Plot of K as a function of t for 0 < t <0.22 and 
log(-K) as a function of t for 0.22 < t < 1.0. = 6 for 

7 = 1.4. K = (%) 2 (1 - M 2 ) = — M = 0, K = 1 , t = 0. 

(l - t )^4 

M = 1, K = 0, t = 1/6; M = 00, K = -00, t = 1. 
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Figure 3.- Plot of K as a function of 
vertical asymptote: a = -0,2513. M 

cr = 0; M = oo, K = -oo, o = -0.2513* 


a. Horizontal asymptote: K = 1; 

= 0, K = 1, a = ooj M = 1, K = 0, 
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Figure 4.- Plot of 
M = 0, 


a as a function of t. 
= 0, cf = coj M = 1, 


a= -0.251251 - (1 - 
= 1/6, a = 0; M = 


1+— + (1_ t) ' 


+ tanh' 


L 3 5 

o, T = 1 . 0 , a = - 0 . 2513 . 
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(a) As functions of t. See equations (l8c) and (22c). M = 0, cd = -oo, 
at t = 0j M = 1; (D = 0 j at t = l/6; M = 1, fl = 0, at t = l/6; 

M = 00 ^ Q = 2.277^ at t = 1.0. 

a) and Q against various functions. 


Figure 5.- Plots of 
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(b) As functions of M. See equations (l8c) and (22c). M = 0, 
M = 1, 0) = 0; M = 1, fi = 0; M = 00 , Q = 2.277. 
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Figure 5.- Concluded 








Figure 6.- Plot of 
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d(log e K 1 / 2 ) 


da) 


as a function of 
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Figure 8.- Plot of € as a function of o> and ft. M = 0, (D = - 00 , 
e = - 00 ; M = 1, co = 0 j € = 0; M = 1, ft = 0, € = 0; M = 00 ^ 
ft = 2.277, € = 2.268. 
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Figure 10.- Different order approximations to curve = — (l - M^)^ 2 

P. 

against o> at free-stream o^,. Horizontal asymptote: kV^ = 1. 


*0 


1/2 = £a(i - 


zero -order approximation by Von Karman; 


= a + bo), first-order approximation; - a + boo + ccd^, 

second-order approximation. 
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Figure 11.- 
in figure 


Other subsonic approximations in 
6. Horizontal asymptote: — = 


comparison with ~-^Log e 
0; vertical asymptote: cd = 0. 
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Figure 12.- Supersonic approximations in comparison with exact 


- 1 against fi. Vertical asymptote: Q = 2.277* 


110 
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Figure 13.- Detailed comparison of third-order approximation with exact 
K(cr). Solid line, exact; triangle, approximation. Horizontal 
asymptote; K = 1; vertical asymptote: a = -0.2513. a = 9.^-, 

b = 1.6l, c = 3 . 89 . % = 

c^ 
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Figure lh.- Comparison of p G /p in hypothetical gas law with exact gas 
law as a function of a. Dotted line, exact; solid line, approximation; 
long dashed line, Loewner approximation (reference 73). Horizontal 

asymptote: ^- = 1. — = O .6339 at a = 0 (M = l). 
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Figure 15. - .Comparison of Mach number of hypothetical gas and exact gas. 
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(a) Physical plane. 


<r 



(b) Hodograph plane (q - 0). (c) Hodograph plane (a - 6). 7 = 1.4. 

a o = 0; = 0.5283; — = 0.6339 

Po P Q 

Figure 16.- Sketch of flow through an inclined-walled straight-edged 

aperture . 
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